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Wc discuss fluctuating order in a quantum disordered phase proximate to a quantum critical point, 
with particular emphasis on fluctuating stripe order. Optimal strategies for extracting information 
concerning such local order from experiments are derived with emphasis on neutron scattering 
and scanning tunneling microscopy. These ideas are tested by application to two model systems - 
the exactly solvable one dimensional electron gas with an impurity, and a weakly-interacting 2D 
electron gas. We extensively review experiments on the cuprate high-temperature superconductors 
which can be analyzed using these strategies. We adduce evidence that stripe correlations are 
widespread in the cuprates. Finally, we compare and contrast the advantages of two limiting 
perspectives on the high-temperature superconductor: weak coupling, in which correlation effects 
are treated as a perturbation on an underlying metallic (although renormalized) Fermi liquid state, 
and strong coupling, in which the magnetism is associated with well deflned localized spins, and 
stripes are viewed as a form of micro-phase separation. We present quantitative indicators that 
the latter view better accounts for the observed stripe phenomena in the cuprates. 
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I. INTRODUCTION 

Ordered states of matter are characterized by broken 
symmetry. Depending on various real-world details, this 
may be relatively easier or harder to detect experimen- 
tally, but once detected it is unambiguous. The notion 
that order parameter fluctuations are important in the 
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disordered phase proximate to an ordered state is a rather 
obvious extension of the idea of broken symmetry; how- 
ever, this notion is more difficult to define precisely. To 
develop this important concept, we will discuss strategies 
for detecting quantum fiuctuating order in the particular 
context of high-temperature superconductor. We will fo- 
cus on the detection of stripe order in the putative stripe 
liquid phase of these systems.^ More generally, we are 
interested in electronic liquid crystalline states and their 
associated fluctuations (Kivelson et al., 1998), but the re- 
sults are also easily generalized to other forms of order. ^ 

Since the whole notion of fiuctuating order is based 
on proximity of an ordered state, it is essential first to 
establish the existence of the "nearby" ordered state by 
directly detecting the relevant broken symmetry. The 
ordered phase may be induced by making small changes 
to the chemical composition of the material, applying 
pressure or magnetic fields, etc. Unless an actual ordered 
state can be reached, it is dangerous to speculate about 
the effects of related fluctuating order. 

Typically, the best way to detect both the broken- 
symmetry state and the relevant fluctuations is by 
measuring the appropriate dynamical structure factor, 
S{q,Lu). Indeed, X-ray and neutron scattering studies 
have provided the best evidence^ of ordered and fluc- 
tuating stripe phases, as we shall discuss. However, 
in many interesting materials, appropriate crystals are 
not available and so such experiments are not possible. 
Here, probes of local order, such as nuclear magnetic res- 
onance (NMR), nuclear quadrupole resonance (NQR), 
muon spin rotation (/^SR), and scanning tunneling mi- 
croscopy (STM) techniques, may be the best available. 
All of these are quasi-static {i.e., nearly zero frequency) 
probes. In a pure quantum system, in its disordered 
phase, the order-parameter fluctuations are not static, 
but rather fluctuate with a characteristic frequency that 
grows with the distance £ to the quantum critical point. 
Thus, unless something is done to "pin" these fluctua- 
tions, they are invisible to local probes. Such pinning 
is induced by boundaries, vortices, crystal-field effects, 
weak quenched disorder, etc. 

In this article, we first discuss results obtained for 
solvable model systems, which we analyze in various ways 
to illustrate the optimal strategies for extracting infor- 
mation about local order. In particular, much of this 
discussion addresses the character of the local order in 
a quantum disordered phase close to a quantum critical 



^ For brief reviews of the current evidence on various types of stripe 
order in the cuprates, see (Carlson et al, 2003; Castro Neto and 
Morais Smith, 200,3: Emery et al, 1999; Orenstein and MiUis, 
2000; Sachdev, 2002; Zaancn, 2000) 

^ Fluctuations associated with off-diagonal {i.e., superconducting) 
order are best studied using different strategies; see, for exam- 
ple, (Carlson et al, 2000; Iguchi et al, 2001; Janko et al, 1999; 
Ussishkin et al, 2002; Wang et al, 2001; Xu et al, 2000) 

^ For reviews of neutron and X-ray evidence of stripes, see (Emery 
et al, 1999; Tranquada, 1998a,b). 



point; we believe that the intuitions gleaned from this 
study are more generally valid, but without the nearness 
to the critical point as a small parameter, it is difficult to 
treat the problem in a controlled fashion. We also review 
in some depth, although by no means exhaustively, the 
experimental evidence of various forms of stripe order in 
the cuprate superconductors. There are several related 
topics that we do not cover in this article; instead, we 
direct the interested reader to recent reviews. Specifi- 
cally, the mechanism of stripe formation^ is given short 
shrift, and the possible relevance^ of local stripe order 
to the mechanism of high temperature superconductiv- 
ity and to the various remarkable normal state properties 
observed in the cuprates is only touched on briefiy. More 
broadly, the context for this review is the role of "com- 
peting orders" in high temperature superconductors.^ 

In Section II we discuss basic scaling considerations 
that govern the relevant length, frequency, and energy 
scales of fluctuating order in the neighborhood of a quan- 
tum critical point, and the consequences of the pinning 
induced by disorder and other perturbations. We show 
that it is the low frequency part of the dynamical struc- 
ture factor, rather than its integral over all frequencies, 
that contains the clearest information concerning local 
order. We also define a new response function which de- 
termines the local density-of-states modulations induced 
by a weak-impurity potential. Finally, we sketch various 
possible versions of the phase diagram of a system with 
competing stripe and superconducting phases. 

In Section HI (and Appendix A), we consider a the- 
oretically well understood system, the one dimensional 
electron gas (IDEG) in the presence of an impurity. We 
compute the local density of states (LDOS), as would be 
measured in STM, and the dynamical structure factor, 
and show by this explicit example how the general con- 
siderations articulated in the present paper are manifest 
in this solvable model. (The IDEG can also be viewed 
as a quantum critical state associated with charge den- 
sity wave (CDW) order.) In addition to its pedagogical 
value, this section contains explicit results that should 
be useful for analyzing STM experiments on quasi ID 
systems with dilute impurities, such as the chain-layers 



* Recent reviews of the mechanisms of stripe formation in the 
cuprates and, more generally in doped antiferromagnets, often 
with rather different perspectives on the issue, can be found 
in Refs. (Carlson et al, 2003; Emery et al, 1999; Hasselmann 
et al, 2002; Ichioka and Machida, 1999; Vojta and Sachdev, 1999; 
White et al, 2002; White and Scalapino, 2000; Zaanen, 1998; Za- 
anen and Littlewood, 1994). 

* Various perspectives concerning the relevance of local stripe or- 
der to other properties of the high-temperature superconductor 
are reviewed in (Carlson et al, 2003; Castro Neto, 2001; Loren- 
zana et al, 2001a,b; Zaanen et al, 2001). 

® Several recent discussions of the role of competing orders in de- 
termining the phase diagram of the high temperature supercon- 
ductors are contained in (Cliakravarty et al, 2001b; Emery et al, 
1999; Kivelson et al, 2001; Orenstein and Millis, 2000; Sachdev, 
2000; Varma, 1997; Zaanen, 2000). 
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of YBCO (Dcrro et at, 2002), or STM experiments on 
carbon nanotubes (Hornbaker et ai, 2002; Odom et ai, 
2002). 

In Section IV we calculate the LDOS in the context 
of a simple model-system with quasiparticles and incip- 
ient order - the weakly-interacting electron gas in two 
dimensions. 

In Section V we discuss applications of these ideas to 
experiments in the cuprates. In Section V.A we dis- 
cuss diffraction studies, and in Section V.B STM studies. 
This latter section contains new insights concerning the 
optimal way to analyze STM data to extract informa- 
tion about fluctuating order; applying these ideas to the 
recent experimental results of Hoffman et al. (2002a,b), 
and Howald et al. (2003, 2002), reveals the existence of a 
nearby stripe ordered phase in optimally doped BSCCO. 
^ Section V.C contains a discussion of experiments to de- 
tect local nematic order, and a discussion of recent STM 
evidence of such order on the same BSCCO surfaces. In 
Section V.D, we discuss indirect evidence of stripe or- 
der that comes from the "1/8 anomaly." In Section V.E 
we briefly discuss other probes of stripe order, including 
NMR, NQR and /iSR. 

In Section VI, we address an important issue of per- 
spective: In the weak coupling limit, the electronic prop- 
erties of a solid are essentially determined by the quasi- 
particle band structure, while collective modes and var- 
ious forms of order reflect slight rearrangements of the 
states near the Fermi surface. Often these effects can 
be plausibly treated in the context of a Hartree-Fock or 
RPA treatment of the residual interactions. Conversely, 
in the strong coupling limit, the physics is more sim- 
ply understood in terms of interacting collective modes, 
such as spin waves, superconducting phase fluctuations, 
and the "phonons" of a charge-ordered state. Here, we 
discuss the interrelation between these opposite perspec- 
tives. In particular we present evidence that the stripe 
order observed in the cuprates and related compounds is 
best thought of as arising from the strong-coupling an- 
tiferromagnetism of the undoped Mott insulating parent 
compounds than from the more conventionally metallic 
physics of the strongly overdoped materials. 

Our most important conclusions are embodied in five 
numbered "Lessons," which are stated in Sections II, III 
and IV, and summarized in Section VII. 



II. GENERAL CONSIDERATIONS 

Since the discovery(Bednorz and Miiller, 1986) of high 
temperature superconductivity in the cuprate supercon- 
ductors, there has been intense interest in the question 



This conclusion is in agreement with the inferences drawn from 
the data by Howald et al. (2003, 2002), but in disagreement with 
those drawn by Hoffman et al. (2002b); we will discuss the origins 
of the differing conclusions. 



of how a Mott insulating antiferromagnet is converted, 
upon doping, into a high temperature superconductor. 
Generally, in electronic systems, there is a competition 
between the kinetic energy (Fermi pressure) which fa- 
vors a uniform Fermi liquid phase with sharply defined, 
itinerant quasiparticles, and the Coulomb repulsion be- 
tween electrons, which favors various forms of insulat- 
ing magnetic and/or charge ordered states. Thus, it 
should not be at all surprising to find various forms 
of charge ordered states appearing in doped antiferro- 
magnets. In particular, "stripe" states have consistently 
turned up in theoretical studies of doped antiferromag- 
nets, including early Hartree-Fock solutions of the Hub- 
bard model^, studies of Coulomb frustrated phase separa- 
tion^, slave-boson mean-field theories of the t-J model^°, 
and various Monte-Carlo and DMRG studies of t-J and 
Hubbard models^^. The discovery of stripe order in 
La2-a;Sra;Ni04-i-5 and soon after in Lai.6-xNdo.4Sra;Cu04 
(Tranquada et al, 1995b) added considerable credibility 
to the suggestion that stripe states form an important 
bridge between the Mott insulator, and the more metal- 
lic state at heavy doping. 

"Stripes" is a term that is used to describe unidirec- 
tional density wave states, which can involve unidirec- 
tional charge modulations ("charge stripes") or coex- 
isting charge and spin density order ("spin stripes" )."'^^ 
CDW states can occur in the weak coupling limit if there 
are sufficiently well nested portions of the Fermi sur- 
face. From the strong coupling perspective, stripes are a 
real-space pattern of micro-phase separation (hence the 
name), driven largely by a lowering of the doped- hole 
kinetic energy in which the doped holes are itinerant 
in metallic rivers, and the antiferromagnetic correlations 
of the parent insulator are preserved in between. Al- 
though the characteristic dynamics of the stripes may 
be quite different in these two limits, as are the possible 
implications for other physical properties, from a bro- 
ken symmetry viewpoint there is no difference between a 
unidirectional CDW and a stripe-ordered state. 

Since the principal purpose of the present paper is 
to address questions concerning the existence of local 



(Kfachida, 1989; Schulz, 1989; Zaanen and Gunnarsson, 1989) 
(Castellani et al, 1995; Emery and Kivelson, 1993; Hasselmann 
et al, 1999; Low et al, 1994; Seul and Andelman, 1995) 
(Han et al, 2001; Lorenzana and Seibold, 2002; Seibold et al, 
1998) 

For varying perspectives, see (Hellberg and ManousaJsis, 1997, 
2000; White and Scalapino, 1998, 2000) 

(Chen et al, 1993; Hayden et al, 1992; Tranquada et al, 1994) 
One can also imagine "orbital stripes" which involve a unidi- 
rectional modulation of a staggered flux or d-density wave or- 
der (SchoUwoeck et al, 2002); in large part, the analysis of the 
present paper would apply equally well to this form of order, as 
well. 

As discussed, for instance, in (Castro Neto and Morais Smith, 
2003; Chernyshev et al, 2000; Emery et al, 1999; Poilblanc and 
Rice, 1989; Schulz, 1989; White and Scalapino, 1998; Zachar, 
2000). 
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stripe order we will for the most part consider unidirec- 
tional CDW and stripe order as two limits of the same 
physics. Indeed, it is important to stress that the micro- 
scopic mechanism of stripe formation in particular ma- 
terials is still not clear. Questions such as whether or 
not the long-range Coulomb interactions are important 
have been widely debated, but remain unresolved. There 
is always a linear coupling between any form of charge 
order and equal period lattice distortions, so it is obvious 
that lattice distortions play a significant role in enhanc- 
ing all the observed stripe phenomena, whether or not 
they are fundamental to the mechanism. This observa- 
tion is further corroborated by the observed stabilization 
of stripes by particular orientational orders of the api- 
cal oxygens in some of the cuprates (Axe and Crawford, 
1994; Biichner et ai, 1994), and by a series of thermal 
transport (Baberski et al., 1998; Sun et ai, 2003) and 
isotope efi'ect (Crawford et al., 1990) anomalies associ- 
ated with the onset of stripe order (Tranquada et al., 
1995b). 



A. Stripe ordered phases 

In 2D, an ordered stripe phase with the stripes running 
in the y direction gives rise to new Bragg peaks in the 
electronic charge scattering (and corresponding peaks in 
the nuclear scattering) at k = ±Qch — (27r/a)(±(5ch, 0) 
and harmonics, where S^h = 1/ ^ch said Xch^i is the charge- 
stripe period. Where spin order coexists with charge or- 
der, new magnetic Bragg peaks occur at harmonics of 
k = Qs = Qaf ± ^Qch where Qaf = (27r/a)(|, |) is the 
Neel ordering vector. 

Charge stripes break rotational symmetry and trans- 
lation symmetry perpendicular to the stripes-in a crys- 
tal, these are to be interpreted as breaking of the crystal 
symmetry group, rather than the continuous symmetries 
of free space. The relevant order parameter, (pq^i,), is 
the Fourier component of the electron charge density at 
the ordering wave vector. If the state is, in addition, 
a conducting or superconducting electron fluid, it is a 
charge stripe smectic. Spin stripes, in addition, break 
spin-rotational and time reversal invariance (although a 
particular combination of time reversal and translation 
is preserved). Wherever there is spin-stripe order, there 
is necessarily (Zachar et ai, 1998) charge-stripe order, as 
well. The new order parameter which distinguishes the 
spin stripe phase, (sqJ, is the Fourier component of the 
spin-density. 

There is more than one possible stripe-liquid phase. In 
particular, there is the possibility of a "stripe nematic" 

phase, in which thermal or quantum fluctuations have 
caused the stripe-ordered state to melt {i.e., translational 



symmetry is restored) but orientational symmetry re- 
mains broken (i.e., a snapshot of the system is more likely 
to see stripe segments oriented in the y than x direction). 
The nematic shares with a charge-ordered state a precise 
deflnition in terms of broken symmetry. The order pa- 
rameter can be taken to be any physical quantity trans- 
forming like a traceless symmetric tensor. For instance, 
the traceless part of the dielectric tensor is a measure of 
nematic order of the charged fluid. In two dimensions, 
a useful nematic order parameter is also (Chaikin and 
Lubensky, 1995) 



5(k) - 5(7^[k]) 



5(k) + 5(7e[k])' 

where 7?, is a rotation by tt/2, and 



(2.1) 



S{^)= f^^SiKu;) (2.2) 

is the thermodynamic (equal-time) structure factor. 
More exotic states can also occur, such as a "nematic 

spin-ncmatic," which breaks rotational, spin rotational, 
and time reversal symmetry, but not the product of TZ 
and time reversal. 



B. Fluctuating order near a quantum critical point 

We consider a system in a quantum disordered phase 
near to a quantum critical point beyond which the ground 
state would be ordered. To be concrete, we will discuss 
charge-stripe order. By the quantum disordered phase, 
we mean one in which there is no long-range stripe order 
as T — » 0, but there may be other forms of order, for 
instance superconducting order. 

The charge-density dynamical structure factor, 
5'ch(k, w), for k in the neighborhood of Qch; measures 
the collective fluctuations which are most sensitive to 
the proximity of the quantum critical point. The scaling 
theory of quantum critical phenomena tells us that, on 
the quantum disordered side of the quantum critical 
point, there arc diverging length and time scales, ^ ^ £~'^ 
and T ^ where is a critical exponent, z is the 

dynamical critical exponent, and £ is the dimensionless 
distance to the quantum critical point. 



Note that, in a crystal, wave vector equalities arc always to be 
interpreted as meaning equal up to a reciprocal lattice vector. 



A nematic spin nematic is most easily pictured as a state in which 
the spin up and spin down electrons each form a nematic state, 
with their principle axes at 90° to each other. 
For extremely cleaj: reviews, see (Sachdev, 2000) and (Sachdev, 
1999a), as well as (Sondhi et al, 1997). 

The values of the critical exponents are "universal" in the sense 
that they are a small discrete set of numbers depending on some 
general features of the critical point. Under many circumstances, 
z is either 1 or 2. i/ depends on the effective dimension, D + z 
and, to a lesser extent, on the symmetry of the order parameter; 
typically, for D > 2, !/ is in the range 2/3 > u > 1/2. 
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In the quantum disordered phase, and in the absence of 
quenched disorder, the fluctuation spectrum has a char- 
acteristic frequency scale, Eq /h. In the absence of dissi- 
pation^^, which is the case we will treat for concreteness, 
Eg typically is a gap in the collective mode spectrum. 
Eg is related to the correlation time by the scaling law 
Eg ~ h/r ~ i"''. For slightly larger than Eg, S'(k, w) 
has a pole corresponding to a sharply defined elemen- 
tary excitation; this "soft-mode" is the quantity which 
condenses across the quantum critical point, so its quan- 
tum numbers and characteristic wave vector directly en- 
code the nature of the nearby ordered state. At some- 
what higher energy (typically, for fkj > 3-Eg) S'(k, w) 
exhibits a multi-particle continuum. Deep in the contin- 
uum the system effectively exhibits scale invariance and 
it behaves in much the same way as if it were precisely at 
the quantum critical point. Both in this high frequency 
regime and at the quantum critical point the notion of a 
quasiparticle is, truly speaking, ill-defined, although for 
systems with small anomalous dimension rj (the typical 
case for D > 2) the continuum will often exhibit fea- 
tures (branch-points, not poles) whose dispersion resem- 
bles that of the Goldstone modes of the ordered phase. 

20 

A classical critical point is described by thermodynam- 
ics alone, so none of these dynamical considerations affect 
the critical phenomena. In particular, it follows from the 
classical limit of the fluctuation-dissipation theorem that 
S'(k) = Tx(k, a; = 0); so 5 has the same critical behavior 
as the (static) susceptibility, Xi even though S involves 
an integral over the dynamical structure factor at all fre- 
quencies. Consequently, a growing peak in S'(k) at the 
ordering vector, k = Qch with width |k— Qch| ~ 1/^, and 
amplitude <S'(Qch) ~ |T — Td"'*' reflects the presence of 
fluctuating stripe order near a thermal transition. Here 
7 = iy{2 — rj) is typically about 1 or greater, so the am- 
plitude is strongly divergent. 

The situation is quite different at a quantum criti- 
cal point, where the dynamics and the thermodynam- 
ics are inextricably linked. Here, the fact that the 
largest contribution to S'(k) comes from the high fre- 
quency multi-particle continuum which does not directly 



In the presence of zero temperature dissipation, such as one ex- 
pects in a normal metal, the collective-mode spectrum may be 
gapless, in which case the interpretation of r is more subtle. 
Most of the scaling arguments we make in the present paper are 
readily generalized to this case, as well. For a recent perspec- 
tive on the effect of dissipation on quantum critical phenomena, 
see (Kapitulnik et al, 2001). See also Millis (1993) and Sachdev 
(1999b) on quantum critical points in metals. 
The naive argument, which is often made, that at frequencies 
larger than Eq the quantum disordered phase looks like the or- 
dered phase is thus not completely correct. There is, however, a 
more nearly correct version of this statement: at high frequen- 
cies, on both the ordered and disordered sides of the quantum 
critical point, the system looks quantum critical, and so looks 
the same whichever state is being probed. 



probe the fluctuating order means that the relevant struc- 
ture in 5(k) is relatively small, and can be difficult 
to detect in practice. To see this, we compare the 
expressions in terms of the dissipative response func- 
tion, X (k, w), for the susceptibility (obtained from the 
Kramers-Kronig relation) and 5(k) (obtained from the 
fluctuation-dissipation theorem): 

X(k,a; = 0)= j^x!'MI^ (2.3) 

^(k)= /£coth(^)x"(k,-) 

In the limit T ^ 0, coth(/3fia>) — > 1; then clearly, the 
factor weights the important low frequency part of 
the integral expression for x much more heavily than 
in the expression for S. Consequently, it follows under 
fairly general circumstances(Chubukov et al, 1994) from 
the scaling form of x" that ^(Qch) ^ '^~"'^x(Qch)j or in 
other words that x(Qch) ~ is strongly divergent in 
the neighborhood of the quantum critical point, while 
S'(Qch) P'"^''' = ^-'^(2-^-')) is much more weakly di- 
vergent, or, if perchance z > 2 — rj, not divergent at 
all. At criticality {£ = 0), similar considerations imply 
that S'(k, cj) is a less singular ftmction at small momen- 
tum/frequency than the corresponding susceptibility. 

These considerations apply even on the ordered side of 
the quantum critical point, where the Bragg peak makes 
an exceedingly small (although sharp) contribution to 
the total structure factor which vanishes in proportion to 
i"^^. By contrast, the Bragg peak constitutes the entire 
contribution to 5(k, to = 0), and is the dominant piece of 
the strticture factor integrated over any small frequency 
window which includes a; = 0. 

Lesson #1: There is an important lesson to be 
gleaned from this general discussion concerning the best 
way to analyze experiments . It is the low frequency part 
of S{k, oj) which is most directly sensitive to the stripe 
order. Rather than analyze S{\s.), information concerning 
local stripe order is best obtained by studying x(k, w = 0) 
or by analyzing the partially-frequency-integrated spec- 
tral function, 

S{K O) = (2fi)-i ^ 5ch(k, a;) (2.4) 
J-n ^Ti" 

where on the ordered side of the critical point, the smaller 
n the better, while on the quantum disordered side, 17 
must be taken larger than Eg /h, but not more than a few 
times Eg / h. This quantity is less severely contaminated 
by a "background" arising from the incoherent high en- 
ergy excitations. As we will discuss in Sec. V.A, this 
is precisely the way local stripe order is best detected in 
clean, high temperature superconductors. 

C. How weak disorder can make life simpler 

In the absence of quenched disorder, there is no useful 
information available from the structure factor at fre- 
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quencics less than Eq, which means in particular that 
static experiments (w = 0) arc blind to fluctuating or- 
der. However, a small amount of quenched disorder, with 
characteristic energy scale Vdis ~ Eq but much smaller 
than the "band width" of the continuum of the spec- 
tral function, has important effects on the low-energy 
states. It is intuitively clear that induced low-frequency 
structure of S'(k, will be largest for values of k where, 
in the absence of disorder, S has spectral weight at the 
lowest frequencies, i.e., for k ~ Q^h- Although the ef- 
fects of quenched randomness on the pure critical theory 
can be subtle (Fisher, 1992; Griffiths, 1969; McCoy and 
Wu, 1968), and are rarely well understood (unless dis- 
order happens to be irrelevant), low-energy states are 
frequently produced. 

Lesson # 2: The upshot is that the effect of weak 
quenched randomness in a quantum disordered phase is 
to produce a low frequency "quasi-elastic" part of the 
spectral function S{k,u>). In other words, disorder will 
eliminate the spectral gap, but will only weakly affect the 
partially integrated spectral function, Eq. (2.4), with the 
integration scale set by Vdis- In particular, in the presence 
of weak disorder, the static structure factor, S'(k, = 0), 
should exhibit similar k dependence as 5'(k, f2) of the 
pure system, and so can be used to reveal the nature of 
the nearby ordered phase. 

1. Response functions 

To formalize some of these notions, we consider the 
somewhat simpler problem of the response of the sys- 
tem to a weak, applied external field which couples to 
the order parameter. In a quantum disordered phase 
the existence of an order parameter can only be made 
apparent directly by coupling the system to a suitable 
symmetry-breaking field. For a charge-stripe smectic, 
a non-uniform potential couples to the order parameter 
and thus serves as a suitable symmetry breaking field. 
Thus, the Fourier component Vk of a weak potential in- 
duces a non-vanishing expectation value of the order pa- 
rameter which, in linear response, is 

(Pk) = Xch(k)li + . . . (2.5) 

The linear response law of Eq. (2.5) is valid provided Vk is 
sufficiently small. However, if the typical magnitude of V 
with Fourier components in the range |k — Qch| ^ ^ 1 is 
not small, V^^h ^ -^'G, then the critical region is accessed 
and Eq. (2.5) is replaced by the law 

iPK.) ~ [K.]'^' (2-6) 

where 6 = 2/{D — 2 + rj) is another critical exponent. 

Many local probes, including STM and NMR, are more 
sensitive to the local density of electronic states, A^(r, E). 
Again, in the quantum disordered state, translation in- 
variance implies that, in the absence of an external per- 
turbation, J\f{r,E) = No(E) is independent of r. In the 



presence of a weak potential, it is possible to define a 
relation for the local density of states 

Af(k,£;) = XDos(k,£;)yk + ... (2.7) 

where A^(k. E) is the Fourier transform of jV(r, E). Prom 
linear response theory it follows that 

XDOs(k,i;) = {2tt)-' J drdtdTe'^'-'^'-'eiT) 

X ([{*t(r,t + r),*,(r,r)},n(0)]) (2.8) 

where 9 is the Heaviside function, 4" J. is the electron cre- 
ation operator, and n = "^t^cr is the electron density 
operator. Note that, despite appearances, E is not a fre- 
quency variable, but is rather the energy at which the 
time independent density of states is measured. A sim- 
ple sum- rule relates Xch to xdos^ 

Xch(k) = J dEf{E)xDOs{KE) (2.9) 

where f{E) is the Fermi function. 

It is also interesting to consider stripe orientational 
(nematic) order. Since Xch(k) is a property of the imi- 
form fluid phase, it respects all the symmetries of the 
underlying crystal; in particular, if rotation by 7r/2 is a 
symmetry of the crystal, then Xch(k) = Xch(''^[k]), so 
no information about incipient nematic order can be ob- 
tained to linear order in the applied field. However, the 
leading non-linear response yields a susceptibility for the 
nematic order parameter, defined in Eq. (2.1), as 

Qk = Jdp xnem(k; p) [Fp - Vnip]][V-p + V_n[p]] 

(2.10) 

This non-linear susceptibility (Xnem contains a four- 
density correlator) divca'gc^s at the nematic to isotropic 
quantum critical point, and we would generally expect it 
to be largest for k « Qch- It is important to note that the 
density modulations, reflecting the proximity of stripe or- 
der, are a first order effect, while the nematic response, 
which differentiates Qch from 7?. [Qch], is second order, 
and so will tend to be weaker, even if the nematic quan- 
tum critical point is nearer at hand than that involving 
the stripe-ordered state. 

D. Phase diagrams 

When there is more than one competing ordered phase, 
the phase diagram can be very complicated, and more- 
over can be quite different depending on various micro- 
scopic details. Nonetheless, it is often useful to have a 
concrete realization of the phase diagram in mind, espe- 
cially when thinking about issues related to order, and or- 
der parameter fluctuations. In this subsection, we sketch 
schematic phase diagrams which characterize some of the 



7 



(a) 



Stripe 
Order , 

/ \ 1 Superconducting ^ 

Flu/tudting Stripes 





at a quantum critical point deep in the superconduct- 
ing phase. In the region marked "fluctuating stripes," 
there is significant local stripe order whose character is 
governed by the proximity to the quantum critical point. 
Notice that this fluctuation region extends into the stripe 
ordered phase itself; although there is true long-range 
stripe order in this region of the phase diagram, close to 
the quantum critical point the degree of striping at the 
local level is much greater than the small ordered compo- 
nent. There is no universal statement possible concern- 
ing how far the fluctuating stripe region extends beyond 
the ordered phase. Clearly, if the region in which there 
is local stripe order does not at least include the entire 
region of the phase diagram in which high temperature 
superconductivity occurs, it cannot be essential for the 
mechanism of superconductivity. 

Under many circumstances, the stripe ordered phase 
will not end at a quantum critical point, but rather will 
terminate at a first order line (Kivelson et al., 2001). If 
the transition is only weakly first order, this will not af- 
fect the general considerations presented here. If it is 
strongly first order, fluctuation effects are much weaker, 
and the effects of even very weak quenched disorder are 
much more severe than in the case of a continuous tran- 
sition. This limit warrants further study (Cardy, 1999), 
both for application to the cuprates and more generally; 
it is likely relevant for the organic superconductors.^^ 

Figure 1 (b) is a more ornate version of the phase dia- 
gram, in which all the different broken-symmetry phases 
discussed in the present paper are exhibited. Here, we 
illustrate the case in which the nematic, charge, and spin 
ordering occur at distinct transitions, under which cir- 
cumstances the nematic ordering must precede the charge 
ordering, which must in turn precede the spin ordering 
(Zachar et ai, 1998). In general, there are also distinc- 
tions between commensurate and incommensurate, diag- 
onal and vertical, bond-centered and site centered stripe 
order, so the same general considerations could lead to 
significantly more complicated phase diagrams. 

Since the superconducting order is suppressed in the 
vortex state, and eliminated above Hc2, studies of the 
magnetic-field-induced changes in the local stripe order 
have recently emerged as one of the best ways of deter- 



FIG. 1 Schematic phase diagrams showing various forms of 
stripe order, and its interactions with high temperature su- 
perconductivity; see text for discussion. 

ordered phases and fluctuation effects discussed above; 
the reader is cautioned, however, that the shape, topol- 
ogy, and even the number of ordered phases may vary for 
material-specific reasons. 

Figure 1(a) represents the essential features we have 
in mind, in the absence of any quenched disorder. Here, 
there is a stripe ordered phase which occurs at low doping 
and low temperature, with a phase boundary which ends 



A first order transition of this sort has been reported in the 
BEDT system (Lefebvre et al, 2000) and in {TMTSF)2PF6 (Yu 
et al, 2002). 

In sketching the global shape of the phase diagram, we have 
included the prejudice that static stripe order, especially static 
spin-stripe order, competes strongly with superconductivity, but 
that a degree of local stripe order is necessary for high temper- 
ature pairing-the latter prejudice is reflected in the vanishing of 
the superconducting Tc as the local stripe order is suppressed 
upon overdoping. However, these prejudices affect only the in- 
terplay between stripe order and superconducting order in the 
phase diagram, and not the central features on which we focus 
in the present article, concerning the character of fluctuating 
stripe order. 
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mining the nature of the interplay between stripe and 
superconducting order. There has been a flurry of re- 
cent papers, both experinientaP"^ and theoretical^^ on 
this subject. Since this subject is fairly involved, and is 
well reviewed in the literature, we will not elaborate here 
on the phase diagram in a field, although we will briefly 
mention some of the salient findings in Sec. V. 

Finally, Fig. 1(c) shows the effect of weak disorder 
on the phase diagram illustrated in Fig. 1(b). Phase 
transitions involving breaking of spatial symmetries are 
generally rounded by quenched disorder. The result- 
ing crossovers should have a glassy character i.e., the 
apparent transition temperature is frequency-dependent. 
Thus, the only true phase transitions in the presence of 
disorder are a superconducting and a spin-glass transi- 
tion, although the spin-glass should have a local stripe 
character - it is (Emery and Kivelson, 1993) a "cluster 
spin glass." Note that quenched disorder generally intro- 
duces some frustration (Zachar, 2000), so where in the 
absence of quenched disorder the spin-freezing tempera- 
ture is large, quenched disorder tends to suppress it. Con- 
versely, where quantum fluctuations have destroyed the 
ordered state, the additional "pinning" effect of quenched 
disorder can lead to a spin-glass phase extending beyond 
the border that the magnetic phase would have in the 
zero disorder limit. 



III. ONE-DIMENSIONAL LUTTINGER LIQUID 

In this section, we show how the general features de- 
scribed above play out in the case of the IDEG. For sim- 
plicity, we present results for a spin-rotation-invariant 
Tomanaga-Luttinger liquid (TLL); in a forthcoming pub- 
lication (Bindloss et al., 2003) we will provide details of 
the derivations, and will also discuss the Luther-Emery 
liquid (LEL), i.e., the spin-gap case. In this section we 
will describe the most salient aspects of the ways fluctu- 
ating order manifests itself in these critical systems. In 
Appendix A we give more details of this theory. (We use 
units in which h= ks = 1.) 

Both the TLL and the LEL can be regarded as quan- 
tum critical systems^''' with dynamic critical exponent 
z = 1. So long as the charge Luttinger parameter is in 
the range < Kc < 1 (which is typically the case for 



Recent experimental studies of induced stripe order in the vortex 
state of high temperature superconductors have been carried out 
in (Hoffman et al, 2002a; Katano et al., 2000; Khaykovich et al., 
2002; Lake et al, 2001, 2002). 

Theoretical studies of the effect of magnetic fields on the com- 
petition between stripe and superconducting order have recently 
been extensively studied; sec, for instance, (Demler et al, 2001; 
Kivelson et al, 2002; Zhang et al, 2002) and references therein. 
For an earlier related discussion, see (Arovas et al, 1997; Zhang, 
1997). 

For reviews of the the theory of the IDEG, see (Emery, 1979; 
Fradkin, 1991; Gogolin et al, 1998; Stone, 1994). 



repulsive interactions), the charge susceptibility of the 
TLL diverges as k ^ 2kp as 

x{2kF + q)^\q\-^. (3.1) 

where fcp is the Fermi wave vector and g = 1 — K^. This 
system can rightly be viewed as a quantum critical CDW 
state. (The CDW fluctuations are still stronger in the 
LEL case, where the susceptibility exponent is replaced 
by — 2.) Thus, it is the perfect laboratory for test- 
ing the validity of the general scaling considerations of 
Section II. 

The charge-density structure factor of the TLL is the 
Fourier transform of the more readily evaluated space 
time structure factor 

/OO pQO 
dt I dx S{x,t)e''"'-''^^ (3.2) 
-OO J — CO 

where the charge correlation function <S(r, t) is given by 

S{r, t) = So (r, t) + [e'^fe^-^sfe (r, t) + c.c] 

+ [e'*'=^''54fep(r,i)+c.c.] + ... (3.3) 

Explicit expressions for Sq, S2k-p, etc. arc given in 
the literature. These expressions can be Fourier trans- 
formed (although in general, this must be done numer- 
ically) to yield expressions for the dynamical structure 
factor. Since the LL is quantum critical, this (and other) 
correlation functions have a scaling form; for example, 
near 2kF {i.e., for small q) 

5(2fcp + ,,a;) = l(^)%,..(^,^) (3.4) 

Vc \VcqJ \Vcq 1 J 

where D is an ultraviolet cutoff, i.e., the band width of 
the TLL. The scaling function $2fcF of Eq. (3.4) depends 

implicitly on the dimensionless parameters Kc and ^ 
with Vc and respectively, the charge and spin velocities. 
(Spin rotation invariance constrains the spin Luttinger 
exponent to be equal to 1.) To exhibit the general 
features we are interested in, while making the Fourier 
transform as simple as possible, wc can consider the limit 
T = and set Vs = Vc = v. Then, for q small 

5(2fcp + ,,.;) ^ 1 f ^) % r^i - 1) \4-^ - ii 

V \vqj yv'^q'^ J \_v'^q'^ 

(3.5) 

where 0{x) is the step function. 

For fixed uo, as a function of q, the dynamic structure 
factor exhibits a multi-particle continuum for —ui/v < 
q < co/v, but it does have singular structure, which can 
be thought as an image of the Goldstone (phason) modes 
one would find were there true CDW order: 

as \q\ Lo/v (3.6) 

The equal-time structure factor can also be readily com- 
puted by integrating the dynamical structure factor; as 
q^O, 

S{2kv +q)^A-A! (lgla/2)^^ (3.7) 
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where A and A' are numbers of order 1 and a is a 
short distance cutoff which we take to be a ~ vp/D. 
As promised, this singularity is much weaker than that 
exhibited by X; and indeed the 2kF component of the 
structure factor remains finite (but not differentiable) as 
q^O. (See Fig. 2.) 
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FIG. 2 The T = charge structure factor, S(k) (solid line), 
and the static susceptibility, x(^i'^ = 0) (dashed line), of the 
TLL with Vs = Vo and Kc = 0.5. 



Wc now turn to the spatial structure in the IDEG 
induced by a single impurity at the origin. For Kc in 
this range, the 2kF component of an impurity potential, 
whose amplitude wc will denote by F, is a relevant per- 
turbation (Kane and Fisher, 1992, 1994) with (bound- 
ary) scaling dimension d = ^{Kc + !)• Thus, there is 
a crossover energy scale Tk oz T'^/('^~^^) such that ex- 
citations with energy large compared to Tk see a weak 
backscattering potential, and are thus only weakly per- 
turbed by the impurity. Conversely, for energies low com- 
pared to Tk the system is controlled (Kane and Fisher, 
1992, 1994) by the fixed point at F ^ oo. However, at 
this fixed point, the high energy cutofi^ is replaced by a 
renormalized cutofi^, D Tk- 

We begin, therefore, by considering the limit F = oo, 
i.e., a semi-infinite system with a; > and the boundary 
condition that no current can flow past a; = 0. For finite 
F, this solution is applicable for all energies \E\ <^ Tk- 
The Fourier transform of the impurity induced LDOS 



N{k,E)= N>{k,E)+N<{k,E) (3.8) 

N<{k,E)= — dt dxg<{x,x;t)e'^^*-'''=^ 
27r J_oo Jo 

can be computed using exact expressions for the appro- 
priate single hole Green function (Eggert, 2000; Eggert 



et al., 1996; Mattsson et al, 1997): 

9<{x,x;t) = 5^(*t(a;,t)^r^(a;,0)) (3.9) 

(7 

^ 5o(x, t) + W^'^-^g^k, (x, t) + e-^2fep.g_2^^ (^^ + . . . 

where go{x,t) is the long wavelength part and g2kF{x,t) 
is the 2kp part; the single electron piece can be com- 
puted similarly from g^, or more simply from the relation 
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FIG. 3 Thermally scaled STM spectra N{k, E) in the neigh- 
borhood of A: = 2A:p - (a) and (c) - and ARPES spectra - (b) 
and (d) - near k = fep. In all cases, Vc/va = 4 and Kg = 1; 
for (a) and (b), = 0.5 (6 = 1/16) while for (c) and (d) 
Kc ~ 0.17 (6 = 1/2). In (a) and (c), the density of states 
oscillations are induced by an impurity scatterer at the origin 
with Tk S> E. In the STM spectra, the real and imaginary 
parts are represented by the solid and dashed lines, respec- 
tively. The curves for E/T = 10 and E/T = 20 are offset by 
1 and 2 respectively, in arbitrary units. 
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between their respective spectral densities, N'^{k,E) = 
e-^/TN>{k,E). We will be interested in the 2fcF com- 
ponent which clearly contains information about CDW 
correlations in this semi-infinite IDEG, and for which 
a general expression in space and time has been given 
by Eggert and others (Eggert, 2000; Eggert et al, 1996; 
Mattsson et al., 1997). As with the structure factor, the 
2kF part of N{k, E) can also be expressed in terms of a 
scaling function $ as (see Appendix A) 

where 6 = (1 — Kc)'^ /8Kc and B is a dimcnsionless con- 
stant; we have left implicit the dependence on the dimcn- 
sionless parameters Kc and Vc/vg. 

The scaling form of N{k,E), given in Eq. (3.10), ex- 
presses the fact that the Luttinger liquid is a quantum 
critical system. The spectrum of the Liittinger liquid 
and the existence of a charge-ordered state induced by 
the impurity dictate entirely the structure of the scaling 
function $(.x,?/): it has a multi-soliton continuum for 
both right and left moving excitations, each with lead- 
ing thresholds associated with one-soliton states carrying 
separate spin and charge quantum numbers and moving 
at their respective velocities, as well as a non-propagating 
feature associated with the pinning of the CDW order by 
the impurity. Thus, the STM spectrum exhibits all the 
striking features of the Luttinger liquid: spin-charge sep- 
aration and quantum criticality, i.e., fluctuating order. 
The beauty of the IDEG is that much of this can be 
worked out explicitly. 

In Figs. 3a and 3c, we show the scaling function com- 
puted numerically for various representative values of the 
parameters. The plots show the real and imaginary parts 
of N{2kp + q, E) at fixed E/T as a function of of the 
scaled momentum, Vcq/T, for \q\ <C k-p- For comparison, 
we also show (Orgad et al., 2001) in Figs. 3b and 3d, 
the single hole spectral function A{kY + q, lo) that would 
be measured in an angle resolved photoemission experi- 
ment (ARPES) on the same system in the absence of an 
impurity. 

There are several things to note about these plots^^: 

1. Because of the quantum critical scaling form, high 
energy and low temperature are equivalent. Note, 
however, that in interpreting the large E/T spectra 
as representative of the low-temperature behavior 
of the system, it is important to remember that the 
thermal scaling of the fc-axis of the figure hides the 
fact that all features of the spectrum are becoming 
sharper as T — > 0; this is made apparent in Fig. 4. 

2. It is clear that there are right-dispersing features 
of the scaling functions characterized by the spin 



Details of the asymptotic analytic behaviors of the singularities 
exhibited in these plots can be found in Appendix A. 
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FIG. 4 Low-temperature {E/T = 100) form of the STM spec- 
tra with k near 2fcF for the same parameters as in Fig. 3; 
a) is for Kc = 0.5 and b) Kc = 0.17. Here, we have ex- 
pressed N{k,E) = \N{k,E)\exp[i(t>{k,E)] where the ampli- 
tude is shown as a soUd Une and the phase as a dashed line. 

and charge velocities. The interference between dis- 
persing features of the ARPES spectrum near kp 
and — fcp can be loosely thought of as giving rise 
to the dispersing features in the STM spectrum; 
indeed, as Kc 1 (the non-interacting limit), 
Re{iV(2fcF + q,E)} cx A{kF - q/2,-E) at fixed 
E > 0. However, it is also clear from the fig- 
ure that the stronger the interactions, the less di- 
rect is the resemblance between N{2kF + q, E) and 
A{kF-q/2,-E). 

3. There is also a very weak feature in the spectrum, 
visible only at quite large E/T (compare Figs. 3 
and 4), which disperses in the opposite direction 
(left) to the main features of the spectra with ve- 
locity —Vc- Because these spectra are shown only 
for fc's in the neighborhood of -|-2fcF (or -|-fcF for 
A{k,E)), there is no symmetry between right and 
left moving excitations. If we showed the spectra on 
a larger scale, there would of course be the mirror 
symmetric spectra near — 2fcF (or — fcp) as required 
by Kramer's theorem. 

4. In the STM spectrum, but not in the ARPES spec- 
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trum, there is a feature near q = which does 
not disperse with increasing energy; this is directly 
related to the pinned CDW order. Note that for 
Kc = 0.5, this feature is weak in Fig. 3, and only 
becomes prominent at very large E/T, as shown in 
Fig. 4; but this is the most important feature of 
the data if one is interested in evidence of pinned 
CDW order. 

Lesson #3: Prom this explicit example we learn that 

dispersing features in an STM measurement that resem- 
ble the interference effects that arise from non-interacting 
quasiparticles do not necessarily imply the existence of 
well defined quasiparticles! 

It is interesting to compare N{k,E) with the LDOS 
averaged over some energy scale, 



N{2kF +q,E) = E- 



-7: 



deN{2kF + q,e). (3.11) 



Note that at T = 0, EN{k, E) (pk) as E ^ oo. How- 
ever, the expression we have used for N{k, E) was de- 
rived for an infinite strength scattering potential, and so 
is only valid up to energies of the order of Tk <C D. This 
integrated quantity is shown in Fig. 5 for representa- 
tive parameters'^ at T = 0, and in the limit as q ^ 0, 
N{2kF + q, Tk) has the asymptotic behavior 



iV(2A;F + g, Tk) oc 



1 

TV 



aTK \ 
Vc J 



2b 



Tk 



s/2 



(3.12) 



Both the integrated induced tunneling density of states 
and the induced tunneling density of states at fixed 
voltage E have a singular behavior of the form q~^^^ 
(See Appendix A). However, the big difference is that 
N{k,E) has dispersing singularities in addition to the 
non-dispersing singularity at fc — > 2fcF, while the singu- 
larity at fc ^ 2kF is the only singular feature of N(k, Tk)- 
Thus, N{k,TK) is more easily analyzed for evidence of 
an almost ordered CDW state. 

At energies E > Tk, the response of the system to the 
presence of the impurity is weak (i.e., proportional to 
r), and can be computed in perturbation theory. For in- 
stance, in the weak-impurity limit Tk the integrated 
response of the system over all energies is dominated by 
the high energy, perturbative regime where at T = we 
find, in agreement with the sum rule of Eq. (2.9), (see 
also Appendix A) 

Ni2kF + q,D)= Xch(2fcF + q)T + ©(F^) (3.13) 



By comparing the weak F behavior of N{2kF - 
(3.13)) with its behavior in the opposite F - 



-q,D) (Eq. 
-> oo limit. 



"^"^ For technical reasons, what is actually plotted in Fig. 5 is 
JV(fc) = limT^.^^/r^dE[l - f{E)]N{k,E), where / is the 

Fermi function, but the distinction between this and N{k,TK) 
is not important here. 
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FIG. 5 Thermally scaled integrated STM spectra N{k) = 

|iV(/c)|e^^''=' ~ N{k,TK) in the neighborhood of A; = 2kY for 
Vc/v^ = 4 and = I with a) = 0.5 and b) = 0.17. 
The solid line is the amplitude and the dashed line is the 
phase. ^® 



Eq. (3.12), we see that since Xch ~ (as q — > 0), 

the integrated LDOS for F small is much more singular as 

q ^ than for large F. This is a necessary consequence 
of the fact that the backscattering impurity potential is 
a relevant perturbation . 

It is also worth noting that the impurity induced LDOS 
at a fixed finite distance x from the impurity 



(3.14) 



is always large at low energies \E\ <^ v/x compared to 
the background DOS of the clean TLL 



Af{E) ~ E 



(3.15) 



The meaning of iV(2fcF + q), the Fourier transform LDOS inte- 
grated up to some high energy scale, depends on this scale. For 
r the high-energy cutoff is the band width D, and N{2k-p+q) 
is the Fourier transform of the induced charge density. Similarly, 
for r — > oo, where the high-energy cutoff is large, N(2kp + q) 
is the Fourier transform of the charge density of the semi-infinite 
system. However, for intermediate values of the backscattering 
amplitude P the partially integrated LDOS N(2k-p + q) cannot 
be interpreted as a modulation of the charge density. 
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This is yet another illustration of the way impurities en- 
hance the low energy effects of fluctuating order. 



IV. TWO-DIMENSIONAL FERMI LIQUID 

In this section, we consider the application of these 
ideas to the case of weakly-interacting electrons in 2D. 
To obtain explicit results, we will consider electrons with 
a quadratic dispersion ek = fi^k^/2m. 

Although, as seen in Eq. 2.8, xdos is a two par- 
ticle correlator, for non-interacting quasiparticles, it 
can be expressed as a convolution (Byers et al., 1993; 
Polkovnikov et al, 2002a) of two single-particle Green 
functions: 

XDOs{KE) = -hmy J^Gik + ci,E)G{ci,E) 

(4.1) 

Therefore, in weakly-interacting systems, it can be an- 
alyzed to obtain information about the single-particle 

spectrum. Explicitly, it is easy to see that for free elec- 
trons in 2D xdos{^,E) is given by 



XoiKE) = 



^fi' (ek - 4S) 



(4.2) 



Here, the subscript is introduced for later convenience 
to signify x in the non-interacting limit, and once again 
9(x) is the Heaviside (step) fimction. For fixed E as a 
function of k, this quantity diverges along curves in k 
space where ek = iE (|k| = 2fcF). Note that the fact 
that xa vanishes for ek < 4i? is a peculiarity of the 2D 
case with infinite quasiparticlc lifetime. In 3D, 



1 

8^ 



2to\ 



3/2 



1 



:ln 



ek 



2VE + ^ 



2VE-^ 



(4.3) 



which also has a singularity as ek 4:E, but is non-zero 
(and positive) for all ii^ > 0. In contrast, as shown in 
Section III, in a ID Luttinger liquid the induced density 
of states has a phase jump as k crosses any of the propa- 
gating or non-propagating singularities. It is straightfor- 
ward to see that both in 2D and 3D a finite lifetime of the 
quasiparticles leads to a rounding of the singularities and 
that in 2D it also leads to a positive induced density of 
states for all E > 0. It is also simple to verify that to non- 
linear order in the external potential, the induced density 
of states in 2D may or may not have a phase jump across 
the singularity depending on the details of the perturbing 
potential. For instance, the induced LDOS produced by 
a single impurity can be computed (from the impurity t- 
matrix) ; it exhibits sign reversal (tt phase shift) across the 
singularity at ek = 4S for a repulsive potential, but not 
for an attractive one. In contrast, as we will show below, 
that any k space structure that arises from proximity to a 
quantum critical point is derived from the susceptibility, 
x(k), which is real and positive, so it always produces a 
signal whose phase is constant. 



Lesson # 4: The density of states modulations induced 
by weak disorder in a non-interacting metal are quite dif- 
ferent in character from those expected from the proxim- 
ity to a CDW quantum critical point, both in that they 
disperse strongly as a function of energy, and the peak 
intensities lie along curves (surfaces in 3D) in k space, as 
opposed to the structure associated with isolated points 
in k space expected from near critical fluctuations. 

The charge susceptibility itself, xo(k), is well known 
from many studies of the 2DEG. It has an extremely 
weak non-analyticity, whenever |q| = 2fcF 



Xo(q) 



m 
2^ 



e{q - 2kF) 



4fc2 



q 



(4.4) 



The inverse Foiirier transform of this meager non- 
analyticity is what gives rise to the famous Friedel os- 
cillations in the neighborhood of an isolated impurity. 



(p(r)) ~ ^ cos(2fcFr). 



(4.5) 



What this means is that, for all intents and purpose, 

the Friedel oscillations are all but invisible in the Fourier 
transform of any conceivable STM experiment on a sim- 
ple metal in D > 1. Some non-trivial method of 
data analysis is necessary instead (Briner et al., 1998; 
Sprunger et al., 1997). 

It is worthwhile considering the effects of interactions 
on this picture. For weak enough interactions, the effects 
can be treated in a Hartree-Fock approximation. Thus, 
if we absorb any interaction-induced changes in the band 
parameters into a renormalized band structure, the only 
change in the above analysis is that the external pertur- 
bation, Vk in Eqs. (2.5) and (2.7) must be replaced by an 
effective potential, Vy ^ Vii + J/kPk, where U (which can 
be weakly k dependent) is the strength of the electron- 
electron repulsion. This leads to the usual RPA expres- 
sion for the susceptibilities of the interacting system, and 
to 



Xch(k) 
Xuosi^E) 



[l-?7kXo(k)]-^Xo(k) 
[l-C/kXo(k)]-^Xo(k,i;). 



(4.6) 



^® In practice, the distinction between an N(k., E) which is peaked 
along curves, indicative of quasiparticle interference effects, and 
an N(k, E) peaJced at isolated "ordering" vectors, k = Q, may 
not always be straightforward to establish in experiment. Con- 
sider the case in which there is an anisotropy of strength a in 
the effective mass of the 2DEG, i.e., = fi^{k^ + aky)/2m. 
In the limit of large anisotropy, a S> 1, when the effects of 
the finite k resolution of actual experiments are taken into ac- 
count, apparent peak-like structures can emerge. If we represent 
the effect of finite resolution by integrating the expression in 
Eq. 4.2 over a range of momenta around different points along 
the ellipse = the integrated expression is peaked near 
k = ^2\J1mE ex, where it is a factor of \/a larger than where 
it is minimal near ^/a\^ = ±2%/2mE iy. 
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(It is the fact that E is a. probe energy, not a frequency, 
which is responsible for the fact that it is simply xa(M-)i 
rather than a frequency dependent factor, which appears 
in the expression for XDOs(k, E).) Not surprisingly, since 
Xo(k) is finite for all k, for small U there is little qualita- 
tive difference between xo(k) and x(k). However, if we 
imagine, as is often done (although we are not aware of 
any reason it is justified) that this RPA expression applies 
qualitatively for larger magnitudes of U, then as a func- 
tion of increasing magnitude of U we would eventually 
satisfy a Stoner criterion for a CDW ?7kXo(k) = 1. Here, 
due to the relatively weak k dependence of xo(k), the 
CDW ordering vector Qch is determined as much by the 
k dependence of J7k as by effects intrinsic to the 2DEG. 

Following this line of analysis, let us consider the be- 
havior of xdos (k, E) in the quantum disordered phase 
close to such a putative Stoner instability, where Xch(k) 
is highly peaked at k = Qch- For k far from the ordering 
vector, both the k and E dependence of Xdos are de- 
termined largely by xo(k, E). However, at fixed voltage, 
XDOs(k, E) will also exhibit a peak at a k « Q^h, with a 
voltage dependent intensity proportional to xo(Qch,-E')- 
Note, particularly, that as long as there is a finite cor- 
relation length associated with the incipient order, i.e., 
so long as the peak in Xch(k) has a finite width, the cor- 
responding peak in Xdos (k, E) will not generally occur 
precisely at Qch- To illustrate this point, here and in 
Fig. 6 we adopt as a simple phenomenological model, 
Xch(k) = A(^A:F)^cxp[-^(|k| - Qch)"^], where ^ is the 
stripe correlation length. From Eq. (4.6), it follows that 
in addition to the singularity inherited from xo (k, E) , 
Xdos(K, E) has a peak at a momentum k which satisfies 



k=Q 



ch 



Ek - 2E 



e^-4E 



Vk log[ek] 



(4.7) 



In short, the peak associated with incipient order is 
weakly dispersing (especially at energies far from any 
quasiparticle resonance condition) but so long as ^ is fi- 
nite, the peak is never "non-dispersing." (In contrast, the 
ID example discussed above is quantum critical, so ^ is 
infinite and the feature associated with fiuctuating order 
is strictly non-dispersive.) 

Although the calculations arc somewhat more involved 
(and therefore must be implemented numerically), the 
same sort of weak coupling analysis can be carried out 
in the superconducting state. Oscillations induced by a 
mean-field with period 4, representing stripes with vari- 
ous internal structures in a d-wavc superconductor, were 
carried out by Podolsky et al. (2002). As expected, the 
stripes induce oscillations in the LDOS with the period 
4, but with energy dependences which refiect both the 
quasiparticle dispersion, and the specific stripe structure 
assumed. Impurity induced oscillations in the LDOS in 
a d-wave superconductor were computed by Byers et al. 
(1993) and Wang and Lee (2003), and the effects of prox- 
imity to a stripe-ordered state, at RPA level, were investi- 
gated by Polkovnikov et al. (2002b). In all cases, interfer- 
ence and stripe-related effects interact in fairly complex 



(k.E) 




FIG. 6 The RPA expression for XDOs(k, £) is plotted as a 
function of momentum at a fixed, low energy E — fj, = .2jj,. A 
small decay rate, V = .025/i, is included to round the singu- 
larity and produce a finite xo{k,E) at all k. We have taken 
the phenomenological form for Xch(k) discussed in the text. 
The figure was plotted with ^ = 5fcF~^ and Qch = 1.8A;f. 



ways that require detailed analysis to disentangle. 
Lesson #5: Depending on what regions of k space 
are probed, the STM spectrum can either be dominated 

largely by band structure effects, or by the incipient 
CDW order. The most singular enhancement of the STM 
signal is expected to occur at energies such that a dis- 
persing feature reflecting the underlying band structure 
passes through the ordering wave- vector. 



V. REGARDING EXPERIMENTS IN THE CUPRATES 

Even where broken symmetry associated with stripe 
order has ultimately been proven to exist, establishing 
this fact has often turned out to be difficult for a number 
of practical reasons. In addition, since quenched disor- 
der is always a relevant perturbation (in the renormal- 
ization group sense), macroscopic manifestations of bro- 
ken spatial symmetries are sharply defined only in the 
zero disorder limit. Nevertheless, the existence of some 
form of order which coexists with superconductivity has 
implications for the phase diagram which can, in prin- 
ciple and sometimes in practice, be tested by macro- 
scopic measurements. A particularly revealing set of 
phenomena occur when the strength of the supercon- 
ducting order is modulated by the application of an 
external magnetic field. Moreover, as discussed in 



For experimental studies, see (Ando, 2002; Ando et at, 1997: 
Hawthorn et al, 2003; Hoffman et al, 2002a; Katano et al, 2000 
Khaylcovich et al, 2002; Lake et al, 2001, 2002; Liu et al, 2003 
Ono et al, 2000; Tyler et al, 1998); for theoretical analysis, see 
(Arovas et al, 1997; Chakravarty et al, 2001b; Demler et al, 
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Sec. V.C, from measurements of macroscopic transport 
anisotropics, electronic nematic order {e.g., point-group 
symmetry breaking) has been identified beyond all rea- 
sonable doubt in quantum Hall systems (Du et al., 1999; 
Lilly et al., 1999a), and very compelling evidence for its 
existence has been reported in the last year in under- 
doped La2-xSra:Cu04 and YBa2Cu306-i-j/ (Ando, 2002; 
Ando et al, 2002b). 

However, most searches for stripe order rely on more 
microscopic measurements, especially clastic neutron 
scattering. One aspect of this that has caused consid- 
erable confusion is that, unless an external perturba- 
tion (such as weak crystalline orthorhombicity) aligns 
the stripes in one direction, one generally finds equal 
numbers of y directed and x directed domains, leading 
to quartets of apparently equivalent Bragg peaks, rather 
than the expected pairs. Fortunately, where sufficiently 
long-ranged order exists, it is possible, by carefully an- 
alyzing the scattering data, to distinguish this situa- 
tion from a situation in which the peaks arise from a 
more symmetric "checkerboard" pattern of translation- 
symmetry breaking. From elastic neutron scattering 
(i.e., measurements of both the magnetic and the nu- 
clear S'(k,aj = 0)), and to a lesser extent from X-ray 
scattering, it has been possible to establish the existence 
of stripe-ordered phases in a wide range of members of 
the lanthanum cuprate family of high-temperature super- 
conductor: Lai.6-xNdo.4Sra;Cu04 over the whole range 
of doped hole concentrations from 0.05 < a; < 0.2, 
La2-xBaxCu04 (Fujita et al., 2002), La2-xSrxCu04 for 
^2 0.02 < X < 0.13, and La2Cu04+5 (Lee et al., 1999) 
(including optimally doped material with a doped hole 
concentration !v 0.15 and a superconducting Tc of 42 K). 
Exciting preliminary evidence of charge-stripe order has 
been recently reported, as well, from elastic neutron 
scattering studies on underdoped YBa2Cu306+j/ (with 
y = 0.35 and = 39 K, (Mook et al., 2002)) and 
optimally doped YBa2Cu306+j; (with y = 0.93 and 
Tc = 93 K, (Mook, 2002)). And, it is worth mention- 
ing, stripe order has also been similarly detected in a 
number of non-superconducting doped antiferromagnets, 
including La2-xSrxNi04+5 and the colossal magneto- 
resistance manganites 

Conversely, it is important to note that so far no evi- 
dence of stripe order, or incipient stripe order, has been 
found in any of the electron-doped cuprate superconduc- 
tors. Indeed, all the low-energy magnetism that has been 



2001; Ivanov et al, 2000; Kivelson et al, 2002; Polkovnikov et al, 
2002b; Zhang, 1997; Zhang et al, 2002). 

Ichikawa et al (2000); Tranquada et al (1997); Wakimoto et al 
(2001) 

Kimura et al (1999); Matsuda et al (2000); Niedermayer et al 
(1998); Suzuki et al (1998); Wakimoto et al (1999, 2001) 
Kajimoto et al (2003); Lec and Chcong (1997); Tranquada 
(1998a); Yoshizawa et al (2000) 

Mori et al (1998a,b); Radaelli et al (1999); Wang et al (2000) 



reported to date (Yamada et al., 2003) is peaked at the 
commensurate ordering wavevector, Qaf, rather than at 
an incommensurate wavevector. This is rather strong ev- 
idence that at least spin-stripe order is absent in these 
materials. This does not rule out a possible role for 
charge inhomogeneity; perhaps the electron-doped ma- 
terials are more prone to form "bubble phases" i.e., crys- 
talline phases with more than one doped electron per unit 
cell (Fogler et al, 1996; Seul and Andelman, 1995). 

We now turn to the core problem: Given a system 
which in the absence of quenched disorder or explicit 
symmetry breaking terms is in an isotropic fluid state, 
how is the existence of substantial local stripe order iden- 
tified? 



A. DifTraction from stripes 

Peaks in S'(k, w) at the characteristic stripe ordering 
vectors indicate a degree of local stripe order. The k 
width of these peaks can be interpreted as an indica- 
tion of the spatial extent of local stripe order, and the 
low frequency cutoff as an indication of the typical stripe 
fiuctuation frequency. So long as there is no spontaneous 
symmetry breaking, ^(k, w) necessarily respects all the 
point-group symmetries of the crystal, and thus will nec- 
essarily always show peaks at quartets of k values, never 
the pairs of k values of a single-domain stripe ordered 
state. 

Low frequency spin fiuctuations with relatively sharp 
peaks at incommensurate wave vectors were detected 
many years ago in inelastic neutron scattering stud- 
ies (Cheong et al, 1991; Thurston et al, 1989, 1992) 
of "optimally-doped" {x « 0.15, ~ 38 K) 

La2-2;SrxCu04. However, not until the discovery (Tran- 
quada et al, 1995b) of "honest" stripe-ordered phases 
in the closely related compound Lai.6-xNdo.4SrxCu04 
was the interpretation of these peaks as being due to 
stripe fiuctuations made unambiguously clear. For in- 
stance, as shown in Fig. 7, the magnetic structure fac- 
tor at low temperature and small but finite frequency, 
hbj = i meV, looks very similar (in absolute magnitude 
and width) in both Lai.6_j;Ndo.4SrxCu04 with x = 0.12, 
where elastic scattering indicating statistically ordered 
stripes has been detected (Tranquada et al, 1995b), and 
in La2-i;SrxCu04 with x = 0.15, where no such static 
order is discernible (Yamada et al, 1995). As discussed 
in the previous section, this is precisely the expected^° 
behavior near a quantum critical point, where presum- 
ably the partial substitution of La by Nd has moved 
the system from slightly on the quantum disordered to 
slightly on the ordered side of a stripe ordering quan- 
tum critical point. "^^ We can now confidently characterize 



This interpretation is further supported by the finding that 
the magnetic ^(k, uj) in near optimally doped La2-a;Sra;Cu04 
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FIG. 7 Comparison of constant-energy scans at hw = 
3 meV through an incommensurate magnetic peak (along 
path shown in inset) for (a) Lai.85Sro.i5Cu04 and (b) 
Lai.48Ndo.4Sro.i2Cu04. Both scans are at T = 40 K > Tc. 
Measurement conditions are described in (Tranquada et al, 
1999b). 



La2-2:Srj:Cu04 over an extremely broad range of doping 
as either being in a stripe-ordered, or a nearly ordered 
stripe-liquid phase. 

An important test of this idea comes from studies of the 
changes in S'(k, to) produced by weak disorder. Specifi- 
cally, in Fig. 8, we compare the low frequency (Matsuda 
et ai, 1993) and elastic (Hirota et ai, 1998) pieces of the 
magnetic structure factor of La2-a;Sr2,Cu04 in the pres- 
ence and absence of a small concentration (1.2%) of Zn 
impurities. (Zn substitutes for Cu.) As one might have 
expected, the Zn slightly broadens the k space structure, 
although not enormously (Kimura et al., 1999; Tran- 
quada et ai, 1999a). Most dramatically, the Zn "pins" 
the stripe fluctuations, in the sense that what appears 
only as finite frequency fluctuation effects in the Zn-free 
material, is pushed to lower frequencies and even to cj = 
by the quenched disorder. 



has scaling properties consistent with its being dominated by a 
nearby stripe ordering quantum critical point; see (Aeppli et ai, 
1997). 

It is certain that one effect of the Zn impurities is to pin the 
charge stripe fluctuations. It may also be that the missing Cu 
spin plays an important role in slowing the spin-fluctuations in 
the neighborhood of the Zn; the potential importance of this 
form of coupling, based on the behavior of a Kondo impurity in 
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FIG. 8 Comparison of magnetic scattering measurements in 
La2-a;Sra;Cu04 with and without Zn; all scans are along 
Q ~ + h, ^,0), measured in reciprocal lattice units, (a) 
Scan at i5 = 1.5 meV and T = 7 K, and (b) difference 
between elastic scans measure at 7 K and 80 K, both for 
Lai.s6Sro.i4Cuo.988Zno.oi204 {n = 19 K) (Hirota, 2000). (c) 
Scan ai E = 2 meV and T = 38 K (Yamada et ai, 1998), 
and (d) elastic scans at T = 1.5 K (circles) and 50 K (trian- 
gles) (Kimura et ai, 1999), for Lai.85Sro.i5Cu04 {Tc = 38 K). 
Measurement conditions were different for each panel; see ref- 
erences for details. 



The issue still remains actively debated whether or 
not various fluctuation effects seen in neutron scattering 
studies of the other widely studied families of high tem- 
perature superconductors, especially YBa2Cu306+y and 
Bi2Sr2CaCu20g+5, can be associated with stripe fluctua- 
tions. We will not review this debate here, but will touch 
on it again in Sees. V.B and V.C, below. 

One question arises (Chen et ai, 2002) concerning 
how we can distinguish a fluctuating stripe phase from 
a fluctuating checkerboard phase (which breaks trans- 
lation symmetry but preserves the symmetry under ex- 
change of X and y). Of course, the strongest indication 
that stripes, rather than checkerboards are responsible 
for the observed fluctuations comes from the presence of 
"nearby" stripe-ordered phases, and the absence of any 
clear evidence of actual ordered checkerboard phases in 
any doped antiferromagnet, to date. However, in theo- 
retical studies(Chen et ai, 2002; Fogler et ai, 1996; Low 
et ai, 1994; Moessner and Chalker, 1996) both types of 
order appear to be close to each other in energy, so this 
argument should not be given undue weight. 

Where both spin and charge peaks are observed, it 
turns out that there is a straightforward way to distin- 
guish. From Landau theory, it follows (Zachar et al., 
1998) that there is a preferred relation between the spin 



a system close to a magnetic quantum critical point, has been 
stressed in (Sachdev et al, 1999). 
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and charge-ordering wave- vectors, Qs + Qg = Qch- For 
stripe order, this means that the spin and charge wave- 
vectors arc paraUel to each otlier, and related (up to 
a reciprocal lattice vector) by the relation 2Qs = Qch- 
However, in the case of checkerboard order, the dom- 
inant spin-ordering wave-vector is not parallel to the 
charge-ordering wave vectors; if Qch = (27r/a)(±(5ch, 0) 
and Q^jj = (27r/a)(0, ±(5ch), the corresponding spin- 
ordering vectors, Qg = Qaf ± ('2'n'/a){^Sch,^Sch) and 
Qs = Qaf ± (27r/a)(5(5ch, — 5<5ch) satisfy the requisite 
identities. Thus, the relative orientation of the spin and 
charge peaks can be used to distinguish fluctuating stripe 
order from fluctuating checkerboard order (Tranquada 
et al, 1999a). 



B. STM measurements and stripes 

STM is a static probe and thus cannot detect any 
structure associated with fluctuating order unless some- 
thing pins it ^"^ . Density or "Friedel oscillations" (Friedel, 
1958) in simple metals produced by the presence of a 
defect are directly related to Fermi-surface-derived non- 
analyticities in the susceptibility, x(k). However, "gen- 
eralized Friedel oscillations" can occur in more diverse 
systems in which the relevant structure in x(k) is not 
directly related to any feature of a Fermi surface. For 
instance, a bosonic superfluid on a lattice close to a sec- 
ond order transition to an insulating, bosonic crystalline 
phase would exhibit generalized Friedel oscillations with 
the characteristic wave length of the roton-minimum - 
these oscillations, in a very direct sense, would image the 
fluctuating crystalline order present in the fluid state. 

There are a few important features of generalized 
Friedel oscillations which follow from general principles. 
If the liquid state is proximate to a highly anisotropic 
state, such as a stripe state, the values of k = Q at which 
X has maxima will reflect the pattern of spatial symme- 
try breaking of the ordered state, but x(k) will respect 
the full point-group symmetry of the crystal unless the 
liquid state spontaneously breaks this symmetry, e.g., is 
a nematic. So, the generalized Friedel oscillations around 
a point impurity in a stripe-liquid phase will inevitably 
form a checkerboard pattern, unless some form of ex- 
ternal symmetry breaking field is applied (Polkovnikov 
et ai, 2002b). 

There is another form of spatial modulation of the den- 
sity of states, one with a period which disperses as a func- 
tion of the probe energy, which is sometimes (incorrectly, 
we believe) referred to in the STM literature as Friedel 
oscillations. This latter eff'ect, which was first demon- 
strated by Crommie et al. (1993), is produced by the 



(HoflFman et al., 2002a; Howald et al., 2003; Hudson et al., 1999; 
Lang et al, 2000; Pan et al, 2000a,b, 2001; Renner et al, 1998; 
Sprunger et al, 1997; Yazdani et al, 1999) 



elastic scattering of quasiparticles of given energy off an 
impurity. The resulting interference between scattered 
waves leads to variations of the local density of states 
at wave vectors Q = k — k', where k and k' are the 
wave- vectors of states with energy E = e(k) = e(k'), as 
determined by the band structure, e(k). Generalized ver- 
sions of these oscillations can occur even when there are 
no well defined quasiparticles, so long as there are some 
elementary excitations of the system with a well-defined 
dispersion relation, as is shown in Sec. HI. 

Thus, in STM studies of cuprates we would expect 
stripe correlations to make an appearance as general- 
ized Friedel oscillations, while quasiparticle-like interfer- 
ence is a distinct phenomenon that could also be present. 
The observation of a checkerboard pattern with a 4a pe- 
riod about vortex cores in Bi2Sr2CaCu208-i-5 by Hoff- 
man et al. (2002a) is provocative evidence for pinned 
charge stripes. These results motivated Howald et al. 
(2003, 2002) to search for similar evidence of stripes in 
Bi2Sr2CaCu208+5 with no apphed field. The discussion 
below is focused on approaches for distinguishing and 
enhancing modulations due to generalized Friedel oscil- 
lations associated with such incipient order. Before con- 
tinuing, though, we should note that compelling evidence 
for a "quasiparticle" interference response has been re- 
ported by Hoffman et al. (2002b) and by McEIroy et al. 
(2003). Those experimental results are not in direct con- 
flict with those of Howald et al. (2003, 2002); however, 
there has been some controversy over interpretation, re- 
garding whether the observed spatial modulations in the 
tunneling conductance are explained entirely by the in- 
terference mechanism, or whether interference modula- 
tions coexist with generalized Friedel oscillations. We 
will return to that controversy in a subsection below. 

Modulations which reflect the spectrum of elementary 
excitations are distinguishable (Polkovnikov et al., 2002a) 
from those related to incipient order in a variety of ways: 
whereas incipient order produces effects peaked near iso- 
lated ordering vectors, Q, the single particle effects are 
peaked along extremal curves in k space which disperse as 
a function of E. Peaks in the Fourier transformed LDOS, 
N{\c,E), produced by incipient order tend to be phase 
coherent while other features either have a random phase, 
or a phase that is strongly energy dependent. Indeed, the 
phase information may be the best way to distinguish the 
consequences of incipient order from interference effects. 
Naturally, in the presence of true long range charge order 
N{k, E) should exhibit sharp (resolution-limited) peaks 
which reflect the charge density modulation in real space. 
This effect has been seen in STM measurements of the 
quasi two-dimensional incommensurate CDW system IT- 
TaS2, lightly doped with Nb, by Dai et al. (1991). 

To demonstrate the effect of the phase, assume we have 
measured 7VL(r, E), the LDOS at a particular bias volt- 
age, V = E/e on a sample of size LxL. Its Fourier trans- 
form, Ni^CkjE), has an arbitrary k-dependent but E- 
independent phase, e"''' ''", which depends on the choice 
of origin of coordinates, Tq. As discussed above, A^L(k, E) 
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can be expected to have contributions from incipient or- 
der (with wave vector Q), and dispersing quasiparticles. 
Integrating the signal over a finite energy window yields 
the quantity N(k,E), in which the contributions from 
incipient order from the entire range of integration add 
constructively, while other features tend to interfere de- 
structively. This mode of analysis is particularly useful 
if an energy window can be found in which none of the 
dispersing features expected on the basis of band struc- 
ture considerations have wave vectors equal to the ex- 
pected ordering wave vector, Q. 

Indeed, Howald et al. (2002) have already demon- 
strated the first half of this point in STM studies of 
a very shghtly overdoped sample of Bi2Sr2CaCu2 08+5 
{Tc = 86). They identified a peak which in N{\<i,E) 
at approximately k — ±Qch ~ ±(27r/a)(0.25, 0) and 
k = ±7^[Qch] = ±(27r/a) (0,0.25) (see Fig. 9), and 
showed that at these wave numbers the phase of A^(k, E) 
is energy independent for E between and 40 meV, at 
which point the amplitude crosses through zero i.e., the 
sign of the signal changes. The constancy of the phase 
implies that the location in real-space of the density of 
states modulation is fixed at all energies, which strongly 
indicates that it reflects pinned incipient order. 

However, as explained above, the existence of incipi- 
ent order can be magnified if we integrate the LDOS as 
a function of energy. In Fig. 9(a) we show data from 
Howald et al. (2003, 2002) in which the values of dl/dV 
were obtained by a Fourier transform of a real space im- 
age at various voltages, V = E/e, on a patch of surface 
of size L = 160 A. The origin of coordinates was chosen 
that A^L(k, E) is real and positive for k— Qch. In the fig- 
ure k is taken to lie along the (1,0) direction and what is 
shown is the real part of Nl (k, E) (the imaginary part if 
generally small and noisy.) In Fig. 9b, we show the same 
data averaged over energy from to E, I /V oc A^L(k, E). 
It is apparent that integration enhances the strength of 
the peak at Qch, and depresses the remaining signal, es- 
pecially when E is smaller than the maximum supercon- 
ducting gap value, Aq ~ 35 meV. Indeed, precisely this 
same integration technique was used previously by Hoff- 
man et al. (2002a) to enhance the magnetic-field induced 
"checkerboard" pattern in vortex cores. 



1. Differing interpretations of the STIVI spectra 

As mentioned in the introduction, there is controversy 
in the literature concerning the interpretation of the STM 
spectra, although there seems to be only relatively mi- 
nor disagreements concerning the data itself. Specifically, 



It is important to emphasize that both the phase and ampHtude 
information encoded in A^^ (k, E) are physically significant. The 
commonly used power spectrum, \N E)\'^ , will not lead to 
destructive interference of propagating features when integrated 
over energies. 
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FIG. 9 The real part of the Fourier transform of the lo- 
cal density of states from an STM study of a near opti- 
mally doped crystal of Bi2Sr2CaCu2 08+i. k is taken along 
the (1,0) direction, in a tetragonal convention, a) shows 
dl/dV oc N{\f.,E) measured at fixed voltage, V — E/e. b) 
shows I /V (X N(k,E), which is the integral of the quantity 
shown in a) from V — to V — E/e. 



Hoffman et al. (2002b) and still more recently McEl- 
roy et al. (2003) have suggested that the peaks seen in 
7V(k, E) can be accounted for entirely in terms of the in- 
terference pattern of sharply defined quasiparticles in a 
d-wave superconductor, with no need to invoke incipient 
stripe order - or any non Fermi liquid character of the 
elementary excitations, for that matter. Let us briefly 
review the line of reasoning that leads to this conclusion. 

As we have seen (Lesson #4), interference between 
quasiparticles in 2D naturally produces ridges in A^(k, E) 
along closed curves in k space, rather than the peaks ob- 
served in experiment. However, as was first recognized 
by Wang and Lee (2003), this result is modified by the 
extreme eccentricity of the contours of constant quasi- 
particle energy in a d-wave superconductor, where the 
dispersion in one direction is vp and in the other is pro- 
portional to Aq. Consequently, with finite experimental 
resolution, at energies less than Aq, one obtains peaks 
in the interference patterns at points in k space which 
connect the tips of the contours that surround each of 
the four distinct nodal points on the Fermi surface. 
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Specifically, if we assume that there are well de- 
fined quasiparticles with energy obtained by solving the 
Bogoliubov-de Gennes equations for a d-wave supercon- 
ductor, then for each energy E there are eight val- 
ues qj (two in the neighborhood of each nodal point) 
which simultaneously satisfy the equations e(qj) — 
(i.e., they lie on the normal state Fermi surface) and 
Aq^. = E(qj) = E {i.e., the quasiparticle creation energy 
is E). Peaks in A^(k, E) will then occur at the various dis- 
tinct values of = — q^ , of which there are 7 distinct 
(up to symmetry) values. For instance, the two wave 
vectors with smallest magnitudes are ki2 = 2qx{l, 0) and 

ki3 = \/2{qx — qy) which come from the inter- 

ference between the state qi = (g^., qy) and, respectively, 
the states at q2 = {—qx,%) and q3 — {qy,qx)- The po- 
sitions of the peaks in Af(k, E) thought to correspond 
to ki2 and ki3 from the work Hoffman et al. (2002b) 
are reproduced in Fig. 10. Given the quasiparticle inter- 
pretation of the STM spectrum, then from the observed 
location of any two distinct peak positions ky at a given 
energy, it is possible to reconstruct the positions of all 
eight locations in the Brillouin zone, q^ which give rise 
to all the expected peaks at this energy. Thus, where 7 
distinct peaks are observed (as are reported by McElroy 
et al. (2003) in some ranges of energy), the quasiparticle 
spectrum is highly over-constrained, and this provides a 
very stringent self-consistency check on the quasiparti- 
cle interpretation. According to McElroy et al. (2003), 
the peaks they observe generally pass this consistency 
check. The experimental case for dispersive features in 
7V(k, E) that are consistent with the quasiparticle inter- 
ference mechanism is quite persuasive. 

The controversy over interpretation revolves around 
the degree to which quasiparticle interference is sufficient 
to explain all of the structure in N{k.,E). Examining 
Fig. 10(a), one can see that for energies below ~ 20 meV, 
ki2 lies close to the expected stripe ordering vector, Qch. 
This is consistent with the experimental measurements 
and analysis by Howald et al. (2002), who had interpreted 
this behavior as evidence for stripe order; in contrast, 
(Hoffman et al., 2002b) argued that the continuity with 
the dispersive signal favors a single mechanism based on 
band structure effects. 

In weighing these alternatives, it is useful to consider 
a second consistency check that the quasiparticle picture 
should satisfy. While, as shown in Eq. (4.1), N(\i,E) is 
actually a measure of a two-particle response function, to 
the extent that it is dominated by single-particle effects it 
should be expressible in terms of a convolution of single- 
particle Green functions, as in Eq. (4.1). Thus, it should 
be consistent with the behavior of the spectral function 
measured in ARPES. In evaluating this connection, sig- 
nificant challenges arise with respect to the sufficiency of 
the quasiparticle interpretation. 

In the first place, ARPES consistently reveals strong 
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FIG. 10 Dispersion of the peaks (colored symbols) in 
(k, E) (from STM measurements on near optimally doped 
Bi2Sr2CaCu2 0g+i ) at wavevectors a) ki (along the (1,0) di- 
rection) and b) k2 (along the (1, 1) direction) from Fig. 4 of 
Hoffman et al. (2002b). The solid lines are from a global fit 
to the ARPES spectrum, as described in the text. 



deviations from a non-interacting line shape. '^^ Well be- 
low Tc, the measured energy distribution curve (EDC) 
at fixed q on the Fermi surface consists of a dispersing 
"quasiparticle- like" peak with a small weight, Z <C 1, and 
a broad, and largely featureless multiparticle continuum 
which contains most of the spectral weight. Moreover, 
the quasiparticle peak is always anomalously broad with 
a distinctly non- Fermi liquid temperature and energy de- 
pendence (Gweon et al., 2001; Orgad et al, 2001). If 
this measured spectral function were used to predict the 
structure of N{k,E), then the quasiparticle interference 
features would be weak (Capriotti et al., 2003) (in pro- 
portion to Z''^), which is something of a surprise given the 
apparent robustness of the experimental features. One 
might also expect the interference features to be broad- 
ened significantly by the short quasiparticle lifetimes at 
the probed energies. Perhaps the most important con- 
sequence is that the intensities of the structures should 
have a dramatic temperature dependence, going to zero 



At very low energies and at low temperatures there is good reason 
to believe, and compelling evidence (although not from ARPES 
or STM) to confirm that the nodal quasiparticles in the super- 
conducting state are well defined and long-lived; see May et al. 
(2000). However, most of the dispersing features seen in STM 
are at higher energies, typically at energies on the order of Ao/2, 
where even at low temperatures there is no evidence of well de- 
fined quasiparticles. 
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as Z goes to zero at T^. This provides a crucial test to 
which we shall return shortly. 

In the second place, Aq determined from the disper- 
sion of the quasiparticle-like peak in ARPES differs sig- 
nificantly, especially at low energies, from that obtained 
by Hoffman et al. (2002b) and McElroy et al. (2003) in 
STM. This is illustrated in Fig. 10, as well. The sohd 
lines in the figure were computed (using the method of 
Wang and Lee (2003)) from a global fit to the measured 
ARPES spectrum, as follows: 1) The location of the 
Fermi surface, e(q) — was computed from the phe- 
nomenological "band structure" that has been obtained 
(Damascelli et al., 2002; Norman et al., 1994) from a fit 
to ARPES experiments. 2) The gap function along the 
Fermi surface was assiimcd to be of the simple d-wave 
form, Aq = ^[cos{qxa) — cos{qya)\, a form that has 
been widely found to fit the measured ARPES spectrum 
in optimally and overdoped Bi2Sr2CaCu208+5, although 
significant deviations from this form are often seen in 
underdoped samples (Mesot et al., 1999). The gap ob- 
served in ARPES experiments corresponds to Aq ^ 35- 
45 meV (Ding et al., 1995; Fedorov et al., 1999; Loeser 
et al, 1997) near optimal doping - for purposes of the 
figure, we have adopted Aq = 40 meV. The solid black 
circle in the figure was obtained from the location of the 
nodal point (which determines the zero energy limit of 
all dispersion curves) as determined from ARPES exper- 
iments; this is a feature of the ARPES spectrum that has 
been looked at in great detail by several groups (Ding 
et al, 1997; Valla et al, 1999). There is agreemenf^o 
with 1% accuracy that the nodal point, where E{(\) =0, 
is q = (27r/a)(0. 195, 0.195). (Note that the phenomeno- 
logical band structure was determined from a fit to older, 
less precise ARPES data, which is why the line in Fig. 
10 does not quite approach the correct zero energy limit; 
manifestly, correcting this discrepancy would only exac- 
erbate the disagreement between STM and ARPES.) 

It is clear from Fig. 10, as originally emphasized 
by Howald et al. (2002), that there are significant dis- 
crepancies between the ARPES and STM results, es- 
pecially below 20 meV. In particular, in the region of 
k space near the expected stripe-ordering wave vector, 
k = Qch ~ (27i'/a)(^, 0), the STM spectrum is consid- 
erably less dispersive than would be expected on the 
basis of the ARPES data; it appears that the peak in 
N{\c,E) largely stops dispersing when it reaches this 
magic wave vector (see for instance the discussion sur- 
rounding Eq.(4.7).)''i 



We thank P.Johnson for pointing this out to us. 
The same discrepancy between ARPES and STM at low ener- 
gies can be seen directly from full microscopic calculations of 
^(k, E) for non-interacting quasiparticles. For instance, A'(k, E) 
was recently computed by Polkovnikov et al. (2002a) for non- 
interacting quasiparticles scattering from a point impurity; it is 
apparent from their Fig. 4 that, at energies of 20 meV and below, 
there is a local minimum in the neighborhood of k = Qch > rather 



It is not clear a priori how much significance one should 
attach to the discrepancies between ARPES data and the 
quasiparticle interpretation of the STM. Both arc highly 
surface sensitive probes, so there is always the issue of 
whether either is telling us anything about the bulk prop- 
erties of the materials, but this worry does not affect the 
comparison between the two sets of measurements. 

It seems significant to us that the most serious quan- 
titative differences between the ARPES dispersion and 
those inferred from a quasiparticle interpretation of the 
STM data occur at low energies, and arc particularly 
pronoiuiced for k's near Qch, where effects of incipient 
stripe order are largest. In this context, it is worth noting 
that in the STM studies of Howald et al (2003) shown 
in Fig. 9, a peak at k = Qch is seen in all the curves 
with < 15 meV all the way to E = 0\ (Hoffman et al 
(2002b) and McElroy et al (2003) do not report spectra 
below about 6 meV. It is possible that, as suggested by 
Howald et al. (2003), that this difference is a consequence 
of the stripe signal being washed out due to the existence 
of many domains in the larger field of view used in the 
experiment of Hoffman et al (2002b).) 

Moreover, in their STM study of the magnetic-field- 
induced structure in vortex-core halos, Hoffman et al. 
(2002a) observed oscillations at precisely these same wave 
vectors, k = Qch. in the energy-integrated LDOS (in a 
12 meV window about voltage) . The amplitude of these 
latter oscillations is very large compared with that of the 
disorder-induced signal, and the dominant contribution 
comes from energies around 7 meV. 

We believe that the issues of interpretation can be 
definitively resolved by studies of the temperature and 
impurity concentration dependence of the signal. Since 
the contribution to A^(k, E) from quasiparticle interfer- 
ence is proportional to Z^, and since Z is observed in 
ARPES experiments to have a strong temperature de- 
pendence and vanish above Tc, it should be very easy 
to quench the interference signal by heating; what is left 
at temperatures approaching Tc cannot be due to quasi- 
particle interference. And, of course, above Tc in over- 
doped samples, and above a suitable pseudogap temper- 
ature in underdoped samples, the superconducting gap 
vanishes, so only the more usual rings can possibly arise 
from quasiparticle interference. Conversely, as discussed 
previously, light Zn doping is known to pin stripes in 
La2-2;Sr2;Cu04, so it is reasonable to expect the same 
effect in Bi2Sr2CaCu208+5. Thus, stripe-related signals 
should be strengthened, and made less sensitive to ther- 
mal depinning by small concentrations of Zn impurities 
in the Cu-0 planes. 

STM studies at elevated temperatures are technically 

demanding, but recently Vershinin et al. (2003) have 
achieved atomic scale resolution at temperatures well 



than a peak! The same is apparent in Figs. 2 and 3 of Wang and 
Lee (2003). 
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above Tc on a Bi2Sr2CaCu208+^ surface. They report 
that for T ~ Tc/2 many of the peaks in A^(k, have, 
indeed, been extinguished. However, the low energy por- 
tion of the peak at ki2 (in Fig. 10) remains unatten- 
uated, and even survives, well above Tc- This effect 
clearly cannot be explained by the quasiparticle inter- 
ference mechanism. An explanation in terms of pinned 
stripes seems much more plausible (although the persis- 
tence of a stripe-induced signal at such an elevated tem- 
perature requires further investigation). 

One further feature of all the data that has been re- 
ported to date is worth commenting on is the large (fac- 
tors of 2 or 3) differences that are observed (Hoffman 
et al., 2002b) in the peak intensities at Qch and i?[Qch]- 
Systematic experiments have not, yet, been conducted 
to verify whether this effect is real, or an experimental 
artifact. If it is not an artifact, then this observation is 
among the first microscopic pieces of evidence of a strong 
local tendency to stripe orientational (i. e., nematic) or- 
der. In addition, such large anisotropics are something 
that cannot be accounted for in any simple way by quasi- 
particle interference (nor local checkerboard order) . It is 
important to bear in mind, however, that even if the ob- 
served anisotropy reflects nematic order, it is expected 
to decrease in magnitude (Howald et al., 2002) in direct 
proportion to A^^^^, where A is the area of the field of 
view, due to the unavoidable domain structure produced 
by quenched disorder. 



C. Detecting nematic order 

While stripe order necessarily implies nematic order, 
the converse is not true. Although nematic order in- 
volves the spontaneous breaking of a spatial (point- 
group) symmetry, when there is no accompanying break- 
ing of translation symmetry, even the identification of 
the ordered state is somewhat subtle, and this holds 
doubly for fluctuation effects. Moreover, since quenched 
disorder is always relevant and results in domain struc- 
ture, true macroscopic measurements of spontaneous ne- 
matic symmetry breaking are not possible. The diffi- 
culty of detecting the rotational-symmetry breaking as- 
sociated with ordered stripes is illustrated by the case of 
La2-xSrxNi04-|-i, a non-superconducting structural ana- 
logue of La2-£cSr2;Cu04. Even transmission electron mi- 
croscopy, which is capable of measuring over a fairly 
small sample area, tends to yield superlatticc diffraction 
peaks for stripe order that reflect the 4-fold symmetry 
of the Ni02 planes in the tetragonal crystal structure 
(Chen et al., 1993). Only recently have electron diffrac- 
tion patterns consistent with the 2-fold symmetry of an 
ordered stripe domain been reported (Li et al, 2003). 
Thus, almost all tests of nematic order in solids neces- 
sarily involve the observation of an unreasonably large, 
and strongly temperature dependent anisotropy in the 
electronic response to a small symmetry breaking field 
which favors the x direction over the y direction. 



1. Transport anisotropies 

An example of a system in which this approach has 

been successfully applied is provided by the 2DEG in 
quantum Hall devices. While the fractional quantum 
Hall effect dominates the physics at very high magnetic 
fields, at intermediate magnetic fields (so that more than 
one Landau level is occupied), it has recently been dis- 
covered (Du et al., 1999; Lilly et al., 1999a) that there 
are a set of anisotropic states which have been identi- 
fied (Cooper et al., 2002; Fradkin and Kivelson, 1999; 
Fradkin et al, 2000) as being quantum Hall nematics, 
and which can be thought of as melted versions of the 
long predicted (Fukuyama et al., 1979; Koulakov et al., 
1996; Moessncr and Chalker, 1996) quantum Hall smec- 
tic (or stripe) phases. These phases are characterized 
by a large resistance anisotropy, p^x ^ Pyy, which on- 
sets very strongly below a characteristic temperature, 
Tn ^ 100 mK. The precise origin of the symmetry break- 
ing field which aligns the nematic domains in these exper- 
iments has not been unambiguously determined (Cooper 
et al, 2001; Zhu et al, 2002). However, by applying 
(Cooper et al, 2002; Lilly et al, 1999b; Pan et al, 2000c) 
an in-plane magnetic field, the magnitude of the symme- 
try breaking field can be varied, and the transition can 
be significantly rounded (giving evidence (Cooper et al, 
2002) that local stripe order persists up to temperatures 
well in excess of Tn), and even the orientation of the ne- 
matic order can be switched, resulting in a state with 
Pxx ^ Pyy 

Experiments that involve such fine control of an ex- 
ternal symmetry-breaking field are considerably harder 
to carry out in the context of the cuprates. Some of 
the relevant materials, such as Bi2Sr2CaCu208+5 and 
La2-xSrxCu04 with x > 0.05, have an orthorhombic axis 
at 45° to the expected stripe directions {i.e., a nematic 
phase can be defined in terms of spontaneous breaking 
of the mirror plane which lies along the orthorhombic 
o-axis.) In these materials, if one wishes to align the 
nematic domains, one must apply a suitable external 
symmetry-breaking field such as a uniaxial strain or in- 
plane magnetic field; however, such experiments tend to 
be challenging. 

When the principle axes of an orthorhombic phase lie 
parallel to the expected stripe directions, the orthorhom- 
bicity (typically < 2% in cuprates) plays the role of 
a small, external symmetry-breaking field. Examples 
where this occurs are superconducting YBa2Cu306+y 
with 0.35 < y < 1 and non-superconducting 
La2-a;Sr2;Cu04 with 0.02 < x < 0.05. In both of these 
cases, resistivity anisotropies as large as a factor of 2 
have been observed (Ando et al, 2002b) in detwinned 
single crystals. Moreover, as in the quantum Hall case, 
this anisotropy is strongly temperature dependent; in 
La2-a;Sra;Cu04, \paa/ Pbb ~ 1| < 10% for temperatures in 
excess of Tn 150 K. Polarization dependent measure- 
ments of the infra-red conductivity on a detwinned x = 
0.03 crystal of La2-a;Sra;Cu04 reveal that the frequency 
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dependence of the conductivity anisotropy has a scale 
comparable to kT^ (Basov, 2002). Large anisotropies 
in the frequency-dependent conductivity have also been 
observed (Basov et al., 1995) in YBa2Cu306+j/ and in 
YBa2Cu408, although in those cases, some part of the 
conductivity anisotropy must be due directly to the Cu- 
O chains. Taken together, these various observations are 
circumstantial evidence of a fair degree of local stripe 
order, and of the presence of a nematic phase; the trou- 
ble being that, without a somewhat quantitative theory 
(which does not exist), it is hard to say how large an 
observed resistance anisotropy must be to be accepted 
as being an "unreasonably large" response to the or- 
thorhombicity, that is, one that can only be understood 
in terms of the alignment of nematic domains. 

More subtle investigations of orientational symmetry 
breaking can be undertaken by studying the transport in 
the presence of a magnetic field. A very clever approach 
along these hues was introduced by Noda et al. (1999). 
They applied a voltage along one axis of the Cu02 planes 
and a magnetic field perpendicular to the planes in order 
to break the 4-fold symmetry of the crystal structure and 
to obtain evidence of one-dimensional charge transport. 
In a stripe-ordered state, the geometry used should be 
sensitive primarily to those domains in which the stripes 
are aligned parallel to the direction of the applied voltage. 
For a; < g, they found that, within the stripe-ordered 
phase, the transverse conductivity tends to zero at low 
temperature while the longitudinal conductivity remains 
finite. (This effect has been explained in terms of the 
electron-hole symmetry of a |;-filled charge stripe (Emery 
et al., 2000; Prelovsek et al., 2001).) 

(Ando et al., 1999) have observed a remarkable 
anisotropy of the resistivity tensor induced by an in-plane 
magnetic field in non-superconducting YBa2Cu306-i-j/ 
with y = 0.32 and y = 0.3. (Presumably, these sam- 
ples are antiferromagiietic.) The results were interpreted 
as evidence for nematic stripe order. Alternative ex- 
planations in terms of anisotropy associated with spin- 
orbit coupling in the antiferromagnetic phase have also 
been proposed (Janossy et al., 2000; Moskvin and Panov, 
2002). This example illustrates the difficulty in making 
a unique association between a bulk anisotropy and mi- 
croscopic stripe order. 

Finally, we note that there is a very direct way to 
detect nematic order using light scattering. This ap- 
proach is well known from studies of classical nematic 
liquid crystals (de Gennes, 1974). Recently Riibhausen 
and coworkers (Riibhausen et al., 2000; Yoon et al., 2000) 
have used light scattering techniques to study the behav- 
ior of the low-frequency dielectric tensor of the mangan- 
ite Bii_xCaxMn03 near and below the charge-ordering 
transition at Tco 160i^. In these experiments, a pro- 
nounced anisotropy of the dielectric tensor was found 
with a sharp temperature-dependence near Tcq. Ra- 
man scattering studies by the same group show that long 
range positional stripe order sets in only at much lower 
temperatures. These experiments suggest that this man- 



ganite is in a nematic state below the transition at Tco, 
and long range charge-stripe order occurs at much lower 
temperatures. 



2. Anisotropic diffraction patterns 

A more microscopic approach is to directly measure 
the nematic order parameter, Qk of Eq. (2.1), by diffrac- 
tion. This was done (more or less) by Mook et al. (2000), 
in neutron scattering studies of the magnetic dynamic 
structure factor of a partially (2 to 1) detwinned sample 
of YBa2Cu306+y with y = 0.6 and « 60 K. No elastic 
scattering corresponding to actual stripe order was de- 
tected; however, well developed structure was observed 
in the inelastic spectrum at the two-dimensional wave 
vectors Qs « (0.5 ± 0.1, 0.5) and « (0.5, 0.5 ± 0.1) in 
units of Itx ja. Remarkably, the intensities of the peaks 
at Qs were found to be about a factor of 2 larger than 
those of the Qg peaks, consistent with the supposition 
that, in a single twin domain, the incommensurate inelas- 
tic structure is entirely associated with ordering vectors 
perpendicular to the chain direction, i.e., Qq k, 1. 

More recently, inelastic neutron scattering studies of 
Stock et al. (2003) on a nearly single-domain sample of 
the Ortho-II phase of YBa2Cu306.5 (Stock et al, 2002) 
have revealed substantial structure down to the lowest 
energies in the magnetic structure factor. These spectra 
are highly anisotropic about the Ncel ordering vector: 
For a scan perpendicular to the chain direction, there is 
a broad flat-topped peak (reminiscent of earlier results 
on underdoped YBCO(Sternlieb et al., 1994)) which is 
strongly suggestive of two barely resolved incommensu- 
rate peaks at the expected stripe ordering wavevector. 
However, a scan along the chain direction reveals a sin- 
gle, sharp peak at the commensurate wave vector 7r/a. 
Moreover, in the normal state, this structure is observed 
at all energies below the "resonant peak" energy, a; w 25 
meV, down to the lowest energies probed. 

Although the presence of chains in this material c;cr- 
tainly means that there is no symmetry operation that 
interchanges the a and h axes, the copper-oxide planes 
are nearly tetragonal. Thus, it seems to us that the ex- 
treme anisotropy of the inelastic scatterings very strong 
evidence of a nematic liquid phase in this material. As 
pointed out by Mook et al. (2000), this conclusion also 
offers a potential explanation for the observed (Basov 
et al., 1995) large but nearly temperature independent 
superfluid anisotropy in the a-b plane. 

Interesting anisotropies have also been observed in op- 
tical phonon branches of YBa2Cu306+j/ with y = 0.6. 
The identification by Mook et al. (2000) of an anomalous 
broadening of a bond-bending mode at a wave vector 
expected for charge-stripe order is potentially the most 
direct evidence of nematic order. However, conflicting 
results have been reported by Pintschovius et al. (2002), 
who have instead observed a zone-boundary softening of 
the bond-stretching mode propagating along b but not 
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for that along a. This latter anisotropy is unlikely to be 
directly associated with a stripe modulation wave vector 
since the anomaly occurs along the direction parallel to 
the Cu-0 chains. Of course, this does not necessarily rule 
out a connection with stripes, as the softening might be 
associated with anisotropic screening due to charge fluc- 
tuations along the stripes. 

In an attempt to understand the effects of stripe or- 
der on phonons, a neutron scattering study (Tranquada 
et ai, 2002) was recently performed on Lai.69Sro.3iNi04. 
Although the charge-ordering wave vector did not play 
an obvious role, the high-energy bond-stretching mode 
propagating parallel (and perpendicular) to the stripe 
modulation exhibited an energy splitting toward the zone 
boundary, while along the Ni-0 bond direction (at 45° to 
the stripes) a softening from zone center to zone bound- 
ary was observed with a magnitude similar to that in the 
cuprates. A better understanding of the nature of the 
relevant electron-phonon coupling processes is required 
to make progress here. 

3. STM imaging of nematic order 

Because it is a local but spatially resolved probe, STM 
is actually the optimal probe of nematic order. One 
way it can be used, which is illustrated in Fig. 11, has 
been explored by Howald et al. (2003, 2002). What is 
shown here is a filtered version of -/V(r, E) measured on 
a patch of surface of a very slightly overdoped crystal of 
BiaSraCaCuzOg+s (T^ = 86 K). Specifically, Howald et 
al. defined a filtered image 

Nf{r,E) ^ J dr'f{r-r')N{r\E) (5.1) 

where in the present case, the filter function, /, has been 
defined so as to accentuate the portions of the signal 
associated with stripe order 

/(r) cx A^e'''^^''^^[cos{Trx/2a) + cos{Try/2a)]. (5.2) 

Clearly, /(r) S{r) when A — > oo, while Nf — 
N{Qch:E) + N{n[Clch],E) + c.c. in the limit A ^ 0. For 
intermediate values of A, the filtered image shows only 
that portion of the signal we have associated with pinned 
stripes. The roughly periodic structure in the image has 
period 4a. We know independently from the analysis in 
Section V.B that there is prominent structure in the raw 
data with this period, but even were there not, the filter- 
ing would build in such structure. However, what is clear 
from the image is that there is a characteristic domain 
structure, within which the stripes appear to dominantly 
lie in one direction or the other. The domain size is 
seen to be on the order of 100 A ~ 25a, which is large 
compared to A~^ and, more importantly, roughly inde- 
pendent (Howald et ai, 2003) of the precise value of A. 
This domain size is a characteristic correlation length of 
the pinned nematic order. 




FIG. 11 A filtered version of the local density of states map, 
Nf(r,E), of the local density of states on a surface patch 
of a Bi2Sr2CaCu20g+i crystal. Here, E — 15 meV and the 
distances on the x and y axes are measured in angstroms. The 
filter is defined in Eq. (5.1) with A = (27r/15a). The arrows 
point along the directions of the Cu-0 bonds. 

This particular method of analysis builds directly on 
the realization of the nematic as a melted stripe-ordered 
state. Indeed, looking at the figure, one can clearly iden- 
tify dislocations and disclinations in what looks like a 
fairly well developed locally ordered stripe array. 

More generally, STM could be used to directly measure 
the two independent components of a suitably defined 
traceless symmetric tensorial density. The simplest such 
quantities are 

QUr,E)^ [dl~dl]N{r,E) 
Q^y{r,E)= 2d.,dyN{r,E) (5.3) 

Of course, these quantities, like the local density of states 
itself, will typically have features which reflect the inter- 
ference between elementary excitations, and other extra- 
neous information. To obtain a better view of the long 
wave-length nematic correlations, we should again inte- 
grate these densities over a suitable energy interval, 
and filter out the short wavelength components: 

^^^""^^ I dr'f{v-v')Q{v',E) (5.4) 

where / might be a Gaussian filter, as in Eq. (5.2), but 
without the cosine factors. A map of Q/ should produce 
a domain structure, similar to that shown in Fig. 11, 
but without all the short wavelength detail. Where the 
domain size, Ljv, is large compared to a, the resulting 
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picture should look qualitatively the same independent of 
the range over which the signal is coarse-grained, so long 
as L^r ^ A~"^ ^ a. The above procedure should work 
well if there are no other long-wavelength features in the 
data. Unfortunately, for BSCCO, the inhomogeneities in 
the gap structure (Howald et al., 2002; Lang et al., 2000) 
hamper such a procedure. 



D. "1/8" anomaly 

Many members of the lanthanum cuprate family of 

high-temperature superconductors exhibit strong singu- 
larities in the doping dependence of various interesting 
low-temperature properties at x = 1/8; together, these 
phenomena are referred to as the "1/8 anomaly." For 
instance, in La2-i:Baj;Cu04 and Lai.6-2;Ndo.4Sr2;Cu04, 
there is (Crawford et at, 1991; Moodenbaugh et al, 
1988) a deep minimum in Tc{x); in La2-xBaa;Cu04 and 
La2-a;Sra;Cu04 there is a pronounced (Crawford et al., 
1990; Zhao et al, 1997) maximum in a(x), the isotope 
exponent a = dlog{Tc)/d\og{M); and in La2-a;Sra;Cu04 
there is a pronounced (Panagopoulos et al., 2002) mini- 
mum in the superfluid density, Us (x) . One of the central 
inferences drawn by Tranquada et al. (1995b) following 
the discovery of stripe order in Lai.g_j;Ndo.4Sr3,Cu04 is 
that this 1/8 anomaly is associated with a commensurate 
lock-in of the stripe structure. At a; = 1/8 the preferred 
spacing between charge stripes is 4 lattice constants, so 
there is an additional commensuration energy which sta- 
bilizes stripe order at this particular hole density. 

While it is possible to imagine other forms of CDW 
order which would similarly be stabilized at x = 1/8, 
there can be no doubt that in this family of materials 
the 1/8 effect is associated with stripe order. This has 
been confirmed, for example, in La2-2;Sra:Cu04 where 
quasi-elastic magnetic scattering from spin-stripe order 
has been detected (Kimura et al., 1999; Suzuki et al., 
1998) for X in the neighborhood of 1/8; in low-energy in- 
elastic measurements, (Yamada et al., 1998) have shown 
that the magnetic peak width is most narrow at a; = 1/8. 
Correspondingly, slow (probably glassy) spin fluctuations 
have been detected by fiSK in the same material with a 
somewhat arbitrarily defined onset temperature, Tg{x), 
which has (Panagopoulos et al, 2002) a pronounced peak 
at X = 1/8. The fact that quasi-elastic magnetic or- 
der as detected by neutron scattering onsets at a con- 
siderably higher temperature than that detected by fiSr 
is clearly a consequence of the inevitable glassiness of 
a density-wave transition in the presence of quenched 
disorder; it reflects the differences in the time scales 
of the two probes, not the presence of two distinct or- 
dering phenomena. More recently, new experiments on 
La2-i;Srj;Cu04 as a function of x show pronounced sin- 
gularities in the x dependence of the c-axis Josephson 
plasma edge (Basov, 2002) and in the low-temperature 
thermal conductivity (Takeya et al., 2002); these effects 
can be interpreted straightforwardly in terms of a peak in 



the stability of the charge-stripe order at x = 1/8. Fur- 
thermore, charge-stripe order has now been detected di- 
rectly by neutron diffraction in Lai.875Bao.i25-a;Sra;Cu04 
(Fujita et al, 2002). 

The large drop in Tc at x = 1/8 found in 
La2-a;Baa:Cu04 is not observed in La2-a;Sr2,Cu04; how- 
ever, such a dip in the doping dependence of Tc can be 
induced in La2-a;Srj:Cu04 (centered at x = 0.115) by 
substitution of 1% Zn for Cu, as shown some time ago by 
Koike et al. (1992) Zn substitution enhances local mag- 
netic order at low temperature near the dip minimum, 
as detected by /iSR (Panagopoulos et al, 2002; Watan- 
abe et al, 2002). Given the clear association between 
the 1/8 anomaly and stripe order in La2-xSrj,Cu04, an 
indirect method to look for the presence of local stripe 
order in other families of cuprate superconductors is to 
test whether a 1/8 anomaly can be induc;cd by Zn substi- 
tution. In the case of YBa2Cu306-i-j/, there already exists 
in the literature strong evidence from the work of Tallon 
and collaborators (Tallon et al, 1995) that the "60-K 
plateau" in the y dependence of is not, primarily, a 
reflection of some form of oxygen ordering in the chain 
layers, as is commonly assumed, h\it is rather a barely 
resolved 1/8 anomaly; this conclusion is also supported 
by a Ca-substitution study (Akoshima and Koike, 1998) 
and by a recent study of transport properties (Ando and 
Segawa, 2002). To further test this idea, Koike and col- 
laborators (Akoshima and Koike, 1998; Akoshima et al, 
2000; Watanabe et al, 2000) studied the doping depen- 
dence of Tc, and of Tg measured by /^SR, in lightly Zn 
doped YBa2Cu;506+j, and Bi2Sr2CaCu208+<5. They have 
found that there is some tendency for light Zn dop- 
ing to produce a dip in Tc(x) and a more impressive 
peak in Tg{x) at x w 1/8 in both YBa2Cu306+j, and 
Bi2Sr2CaCu208+5. Related results have been reported 
(Ando et al, 2002a; Yang et al, 2000) for the single-layer 
system Bi2Sr2_a;Laa;Cu06+5. While it is clear that more 
work is needed to test the connection between the "1/8 
anomaly" and stripe pinning in these various systems, we 
consider this a promising approach to the problem, as it 
permits evidence of local order to be obtained using a 
variety of probes which can be applied in materials for 
which large crystals, and/or easily cleaved surfaces are 
not easily obtained. 



E. Other probes 

It is clear that the most direct evidence for stripes 
comes from techniques that can provide images of 
charge, spin, and/or lattice modulations in real space 
(STM) or in reciprocal space (diffraction techniques). 
We have already discussed neutron and X-ray diffrac- 
tion, but explicit mention should also be made of 
transmission electron microscopy (TEM). The charge 
stripes in La2_a;Sra;Ni04+5 were first detected by TEM 
(Chen et al, 1993), and a recent study has provided 
high-resolution TEM images of local stripe order in 
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Lai.725Sro.275Ni04 (Li et at, 2003). So far, TEM stud- 
ies have not provided positive evidence for stripes in any 
cuprates, but, though challenging, it should be possible 
to do so. 

Less direct but extremely valuable information comes 
from techniques that are sensitive to local order. We 
have already mentioned evidence for local, static hypcr- 
fine fields and slowly fluctuating hypcrfinc fields obtained 
by /xSR. Besides providing a practical measure of local 
magnetic order as a function of doping (Klauss et al., 
2000; Niedermayer et al, 1998), ^SR can detect the dis- 
tribution of local hypcrfinc fields in a sample with rela- 
tively uniform order (Nachumi et al., 1998), as well as be- 
ing able to detect inhomogeneous magnetic or supercon- 
ducting order (Savici et al, 2002). Related information 
can be obtained by NMR and NQR techniques, where 
the latter also provides information on the distribution 
of local cJcK'tric-ficld gradients. Some of the first evidence 
for spatial inhomogeneity of magnetism in lightly-doped 
La2-2:Sr2:Cu04, consistent with stripe-like behavior, was 
obtained in a La NQR experiment by Cho et al. (1992). 
Later experiments (Goto et al., 1997, 1994; Ohsugi et al, 
1994; Tou et al, 1993) (including La and Cu NMR and 
NQR) provided evidence for local magnetic order near 
X = 1/8 in La2-a;Sra;Cu04 and La2-xBaa;Cu04. 

The work by Imai's group (Hunt et at, 1999) suggest- 
ing that local stripe order in variants of La2-a;Sra;Cu04 
could be detected through the "wipeout effect" of Cu 
NQR has motivated a considerable number of further 
studies using Cu and La NMR and NQR (Curro et al., 
2000a; Hunt et al, 2001; Julien et al, 1999, 2001; Simovic 
et al, 2002; Singer et al, 1999; Suh et al, 1999, 2000; 
Teitel'baum et al, 2001, 2000). Although there has been 
some controversy over details of the interpretation (and 
in certain cases, there may be differences between sam- 
ples) , it is now generally agreed that NQR and NMR are 
sensitive to the onset of slow spin fluctuations whose ap- 
pearance tends to correlate with the onset of local charge 
stripe pinning as determined by diffraction. In particular, 
the glassy nature of the ordering has been investigated. 
In the case of Eu-doped La2-a;Sra,Cu04, related infor- 
mation has been obtained by electron-spin resonance de- 
tected from a very low density of Gd impurities (Kataev 
et al, 1997, 1998). A recent Cu NQR and Y NMR study 
(Singer and Imai, 2002) has shown that similar signatures 
are observed in Ca-doped YBa2Cu306, consistent with 
the fiSR study of (Niedermayer et al., 1998), and sug- 
gestive of the presence of pinned stripes in that system 
at low temperatures. Direct evidence for local spatial in- 
homogeneities in La2-a;Sra;Cu04 has been obtained from 
studies of NMR line broadening (Haase et al., 2000) and 
frequency-dependent NQR relaxation rates (Singer et al., 

2002) . 

Very recently, a fascinating study (Haase and Slichter, 

2003) of the NMR/NQR spectra in YBa2Cu306+j, has 
found evidence of two distinct planar O environments in 
the unit cell, but only one Cu environment. This finding 
is not consistent with any form of translation-symmetry 



breaking. It is, however, suggestive of substantial ne- 
matic order, in which the O sites midway between two 
copper sites in one direction (for instance, the chain di- 
rection) has substantially larger hole-density than the O's 
on the bonds in the perpendicular direction. Of course, 
because the crystal is orthorhombic, this does not truly 
imply any symmetry breaking, but it should be possi- 
ble, in principle, to establish whether the magnitude of 
the effect is out of proportion with the small distortions 
produced by the orthorhombicity. 

Several recent NMR studies (Curro et al., 2000b; 
Kakuyanagi et al, 2002a; Mitrovic et al, 2001, 2002) 
have exploited the magnetic-field dependence of the tech- 
nique to probe the spatial variation of nuclear spin-lattice 
relaxation rates in the vortex lattice state of cuprate su- 
perconductors. Vortex cores are regions of suppressed 
superconductivity, so if there is a competing spin-stripe 
ordered phase, it will be enhanced and pinned (Demler 
et al, 2001; Kivelson et al, 2002; Zhang et al, 2002) 
in the neighborhood of the vortex core. Indeed, ^^O 
NMR measurements on YBa2Cu306-i-j/ and YBa2Cu408 
indicate antiferromagnetic-like spin fluctuations and a re- 
duced density of states associated with the vortex cores 
(Kakuyanagi et al, 2002a; Mitrovic et al, 2002). (Sim- 
ilar results, along with evidence of static spin ordering 
in the vortex cores at low temperatures, have been re- 
ported (Kakuyanagi et al, 2002b) in Tl2Ba2Cu06+5.) 
A /iSR study of the details of the magnetic-field dis- 
tribution in YBa2Cu306.5o suggests the presence of a 
small local hyperfine field near vortex cores (Miller et al, 
2002). For completeness, we note that while one neutron 
scattering study has indicated very weak, field-induced, 
elastic antiferromagnetic scattering in superconducting 
YBa2Cu3 064-j, (Vaknin et al, 2000), two other studies, 
focusing on the effect of a magnetic field on inelastic scat- 
tering, found no low frequency enhancement (Bourges 
et al, 1997a; Dai et al, 2000). 

It has been demonstrated that stripe order has an 
impact on phonon heat transport (Hess et al, 1999). 
The thermal conductivity increases slightly on cooling 
through the charge-ordering transition, exhibiting a nor- 
mal peak at ~ 25 K. The suppression of the latter 
peak in superconducting La2-2;Sr2;Cu04 has been at- 
tributed to the scattering of phonons by fiuctuating 
stripes (Baberski et al, 1998). The doping dependence 
of the thermal conductivity measured in YBa2Cu306-i-j/ 
and in mercury-cuprates has been interpreted as evi- 
dence for a 1/8 anomaly in those materials (Cohn et al, 
1999). An enhancement of the thermal conductivity in 
Lai.88Bao.i2Cu04 was originally noted by Sera et al 
(1990) 

Indirect evidence of some form of local charge order can 
also be gleaned from the response of a system to electro- 
magnetic radiation.*^ For instance, Raman scattering 



For a discussion of the electromagnetic signatures of CDW order, 
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(Blumberg et ai, 2002) was recently used to reveal charge 
ordering, and more importantly (for present purposes) to 
detect the effect of collective CDW motion at higher tem- 
peratures in Sri2Cu2404i, a ladder system with a local 
electronic structure very similar to that of the cuprate su- 
perconductors. Careful measurements of the optical re- 
sponse of optimally doped YBa2Cu306-i-j/ have revealed 
(Bernhard et al., 2001; Homes et al, 2000) infrared active 
phonons with oscillator strength comparable to those in 
the undoped insulating compound. This has been plausi- 
bly interpreted as giving evidence of "fluctuating charge 
inhomogeneities" in the copper-oxide planes - the point 
being that for the phonons to be unscreened at that fre- 
quency, their local environment on some scale must be 
insulating. 

ARPES experiments can also be interpreted as giving 
indirect evidence of local stripe order. At an empirical 
level, this can be done by looking (Zhou et ai, 1999, 2001) 
at the evolution of the ARPES spectra in a sequence of 
materials in the La2Cu04 family in which various types 
of stripe order have been detected; this provides a basis 
for identifying similar features in the ARPES spectrum 
of materials in which stripe order has not been estab- 
lished. For instance, stripe long-range order produces a 
tendency for the Fermi surface in the antinodal region 
of the Brillouin zone to be flat (one-dimensional), and 
to suppress spectral weight in the nodal region. Some 
of these same features, especially the differing impact 
on the nodal and antinodal regions, are induced by Zn 
doping of Bi2Sr2CaCu208+5, as well, plausibly reflect- 
ing (White et al, 1999) the tendency of Zn to pin local 
stripe order. More broadly, many of the most striking 
features of the observed (Campuzano et al, 2003; Dam- 
ascefli et al, 2002; Fedorov et al, 1999; Ino et ai, 1999; 
Valla et al., 1999) ARPES spectra can be naturally inter- 
preted in terms of an underlying quasi one-dimensional 
electronic structure (Carlson et ai, 2003). Examples of 
this are the disappearance of the the coherent quasiparti- 
cle peak (Carlson et al., 2001) upon heating above Tc, the 
dichotomy between the widths of the peaks as a function 
of momentum and energy in the normal state spectrum 
(Gweon et ai, 2001; Orgad et al., 2001), and the general 
structure of the spectra and the evolution of the Fermi 
surface with doping (Eroles et al., 2001; Granath et al., 
2001, 2002; Ichioka and Machida, 1999; Salkola et al., 
1996; Zachar et ai, 2000, 2001). 

Finally, several probes have been used to search for 
the inhomogeneous distribution of bond- lengths expected 
from the lattice response to local charge stripe order. 
Pair-distribution-function (PDF) analysis of scattering 
data and extended X-ray absorption fine structure (EX- 
AFS) spectroscopy can detect instantaneous distribu- 
tions of nearest-neighbor bond lengths. Both techniques 
have the potential to provide indirect evidence of dy- 
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namical, as well as static, stripes; however, whether 
they have sufficient sensitivity to detect the very small 
bond-length variations associated with charge modu- 
lation in the cuprates remains a topic of some con- 
tention. (Bozin et al., 2000) reported a doping-dependent 
broadening of the in-plane Cu-0 bond-length distribu- 
tion in La2-xSra;Cu04 detected by PDF analysis of neu- 
tron powder diffraction data. While the reported dop- 
ing dependence is intriguing, the maximum enhance- 
ment of the bond-length spread is much too large to 
be compatible with the Debye- Waller factors measured 
in a single-crystal neutron diffraction experiment on 
Lai.85Sro.i5Cu04 by Braden et al. (2001) A much more 
extreme discrepancy occurs with the EXAFS studies of 
La2-a;Sra;Cu04 (and other cuprates) by Bianconi's group 
(Bianconi et al., 1996; Lanzara et al., 1996; Saini et al., 

1997) . They report a low-temperature splitting of the 
in-plane Cu-0 bond lengths of 0.08 A, which is not only 
incompatible with a large number of diffraction studies, 
but is also inconsistent with EXAFS analysis performed 
by other groups (Haskel et al., 2000; NiemoUer et al., 

1998) . Large atomic mean-square displacements detected 
in YBa2Cu306+j, by ion channeling (Sharma et at, 2000) 
have been attributed to transitions associated with dy- 
namic stripes, but the magnitude of the displacements 
and the temperature dependence are difficult to recon- 
cile with other experiments. 

In short, these various local measures of the distribu- 
tion of lattice displacements are, in principle, a very good 
way to look for local charge-ordering tendencies, but it is 
important to reconcile the results obtained with different 
techniques. Where apparent inconsistencies exist, they 
must be resolved before unambiguous inferences can be 
made. 



VI. WEAK AND STRONG COUPLING PERSPECTIVES 

Order in simple metals is typically thought of as aris- 
ing from a Fermi surface instability caused by the weak 
residual interactions between the lowest energy quasipar- 
ticles. Specifically, density wave order occurs when the 
band structure is such that there arc nearly nested seg- 
ments of the Fermi surface. Thus, even at temperatures 
above the ordering temperature, or if the interactions 
are slightly weaker than the critical strength needed for 
ordering, there will often be structure in the appropri- 
ate dynamical structure factors due to this near nesting 
which in some sense refiect the fact that the metal is close 
to an ordered state. In this section, we will analyze the 
connection between this weak-coupling Fermi-liquid the- 
ory based perspective and the description in terms of the 
collective modes of the ordered or nearly ordered state 
which we have taken until now. 

It is worth mentioning at the outset that one difficulty 
arises from the relative paucity of established theoretical 
results in the strong coupling limit. More than half a 
century of concerted work has produced a rather com- 
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plctc understanding of the effects of weak residual inter- 
actions on the properties of a well formed Fermi liquid in 
a simple metal. However, only in special circumstances is 
a comparably sound and complete understanding avail- 
able in a strongly interacting system. For instance, in 
the weak coupling limit, considerable useful information 
about the quasiparticle spectrum can be adduced directly 
from transport data, because Boltzmann transport the- 
ory relates the two properties in a precise manner. The 
theory of transport in a strongly interacting system is, in 
our opinion, not established. The upshot of this is that 
it is much easier to point to experiments that show that 
simple, weak coupling ideas cannot be safely applied in 
many highly correlated materials, than to point to exper- 
iments that show that strong coupling notions can be. 



A. Distinctions in principle 

1. Thermodynamic distinctions 

As mentioned in the introduction, stripe-ordered states 

are ultimately defined by a specific set of broken sym- 
metries. From this viewpoint, there is no distinction be- 
tween the strong coupling limit, in which spin stripes can 
be viewed as micro phase separation of charges into an 
array of rivers with strips of nearly insulating antiferro- 
magnet between them, and the weak-coupling Hartree- 
Fock description of a SDW which opens small gaps on 
various nearly nested segments of the Fermi surface. In- 
deed, in many circumstances, the two limits are adiabat- 
ically connected as a function of the interaction strength. 
However, this is not guaranteed. For instance, a stripe- 
ordered state can be cither metallic or insulating. At 
T = 0, metallic and insulating stripes are thermodynam- 
ically distinct states of matter. 

There is a more subtle possible distinction possible be- 
tween metallic spin-stripe states in the weak and strong 
coupling limits. Certainly, in the weak-coupling limit, 
the ungapped portions of the Fermi surface support well- 
defined quasiparticles, so the system remains a Fermi 
liquid unless, at still lower temperatures, it suffers an 
additional ordering transition, say to a superconducting 
state. On the other hand, in the strong coupling limit, 
the metallic behavior in the stripe ordered phase may 
or may not be well characterized by Fermi liquid power 
laws. It is now established that non Fermi liquid states 
can, in principle, exist in more than one dimension.*^ 

The character of charge-stripe order tends to more eas- 
ily distinguish between strong and weak coupling limits. 
For repulsive interactions, weak coupling Hartree-Fock 
theories always produce charge order which is parasitic 



In particular the existence of a non Fermi liquid "sliding" phase 
in a stripe ordered system has been shown(Emery et ai, 2000; 
Vishwanath and Ccirpentier, 2001) for a class of interacting Lut- 
tinger liquids. 



on the fundamental spin order. In Landau theory, there 
is a cubic coupling between CDW and SDW order which, 
whenever the primary order parameter is the SDW order, 
leads to CDW order which onsets at the same Tg, but 
with a strength proportional to the square of the SDW 
order: (p2q) oc |Sq|2 ~ |Tc - Tp/^. By contrast, in a 
strong coupling picture of micro-phase separation, if the 
spin-ordering is only triggered when the holes agglom- 
erate into stripes, one would typically expect the stripe 
ordering transition to be first order, or for the charge 
ordering to precede spin ordering. Thus, the measured 
sequence of transitions provide sharp criteria that can be 
used to discriminate between weak and strong coupling. 

In the thermal or quantum disordered phase proximate 
to a stripe ordered state, sharp thermodynamic distinc- 
tions between the strong and weak coupling limits are 
harder to draw. The only possible thermodynamic dis- 
tinction is, again, that in the weak-coupling limit the 
disordered state is necessarily a Fermi liquid, while in 
strong coupling it may either be a Fermi liquid or a non 
Fermi liquid. 



2. Quantitative distinctions 

Even where no sharp thermodynamic distinctions ex- 
ist, there are very many clear physical differences between 
the weak and strong coupling limits. Of these, the most 
obvious is the relative importance of electron quasipar- 
ticles with a sharply defined Fermi surface. In the weak 
coupling limit, the quasiparticle is the essential building 
block in terms of which all other properties are derived. 
In the strong coupling limit, by contrast, even where it 
might be the case that well-defined quasiparticles (with 
very small weight, Z) are recovered in the strict T ^ 
limit, such coherence may require such low temperatures 
as to be only of academic interest. Thus, on a practical 
level, "effective non-Fermi liquids," in which the major- 
ity of the electronic excitations in the experimentally rele- 
vant ranges of temperature, energy, and wave-number arc 
not well defined quasiparticles, cannot be sensibly treated 
in terms of a weak coupling approach. Conversely, where 
well defined quasiparticles are everywhere manifest, there 
is an a priori reason to prefer a weak coupling viewpoint. 

Another set of important of quantitative distinctions 
concern the k and lo dependences of the magnetic dy- 
namical structure factor, S'(k, cj). In a one band- model 
of non-interacting electrons, the integral of S over all k 
and u) is (1/8)[1 — a;^] while in the large interaction (local 
moment) limit it is (1/4)(1 — x), which is not dramati- 
cally different. However, in weak coupling, the integrated 
intensity comes from a very broad range of energies, of 
order Ep, and momenta, of order 2k-F, while for an an- 
tifcrromagnct the integral is dominated by energies less 
than the exchange energy, J, and is highly peaked near 
the magnetic ordering vector, whether or not the system 
is actually ordered. Put another way, in weak coupling 
the structure factor is dominated by the particle hole con- 
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tinuum, and there are signatures of well defined collective 
modes (spin waves) only at energies w < T^,, while in the 
strong coupling limit, deep in the ordered phase, spin- 
waves dominate the structure factor even at energies of 
order J. Also, in the ordered state at weak coupling, the 
ordered moment m at T = is small, m ^ ^b{Tc/ Ep), 
while in the strong coupling limit, not too close to a quan- 
tum critical point, m ^ fiB is large. 

As a concrete example, let us first consider the struc- 
ture factor of a non-interacting electron gas in 2D. For 
momenta g' > 2A;f it is 



k'jp 



2kF vpq 



(6.1) 



when vpqi^^ — 1) < lu < vpqi^^ + 1) and zero oth- 
erwise. (The expression for q < 2kF is somewhat more 
complicated.) Clearly, the spectral weight is distributed 
smoothly over the entire support of this particle-hole 
continuum, i.e., over a range of momenta of order 2A;f 
and of energies of order the band width. From the weak 
coupling perspective the chief qualitative effect of inter- 
actions is to produce (quasi)bound states that manifest 
themselves as peaks at or below threshold of the contin- 
uum. In the broken- symmetry phase these bound states 
are Goldstonc modes, which descend down to zero energy 
and loose weight to the new Bragg peaks. What is im- 
portant to note here is that even in the broken-symmetry 
phase, most of the spectral weight remains in the inco- 
herent background spread more or less uniformly over 
energies of order band width. 

In contrast, consider the zero temperature structure 
factor of an insulating antifcrromagnct computed to first 
order in a 1/S expansion, whore S is formally the magni- 
tude of the spin, but more physically should be viewed as 
the "distance" of the system from the nearest quantum 
disordered phase. Here, on a square lattice, 

Sij (q, w) = S{S - A) rhi m, 5{io) (5(q - Qaf) (6.2) 



5* /2 - cosg^, - cosgy - - \ u ^ 

-\- — \ - [Oij — mimj) oiuj — uja) 



where = 2<S'J-\/4 — (cos kx + cos kyY is the spin wave 
dispersion (which depends on the lattice and the details 
of the exchange interactions), m is a xvaii vector parallel 
to the magnetization, and A = 0.3932+. . .. (Higher order 
terms in 1 /S typically produce finite lifetimes for all but 
the lowest energy spin waves, and also give rise to a multi 
spin wave continuum.) As promised, the basic energy 
scale is sot by J, and a largo portion of the total spectral 
weight is found in a small region of k space around k = tt. 

Note that the energy scale J can be deduced from the 
value of the spin wave energy near the zone boundary, 
Wko = 45 J = 2 J, for ko = (tt, 0). Even in the thermally 
disordered state, so long as T <C J, the zone edge spin 
wave can be studied to get an idea of the relevant energy 
scales for magnetic excitations. 



Finally, distinctions can be drawn based on the sen- 
sitivity of the various signatures of order or proximate 
order to various weak perturbations: weak disorder, as 
we have discussed extensively above, can serve to pin 
fluctuating order, and hence enhance the small w por- 
tion of 5(q, Lo). However, in weak coupling, the principal 
effect of weak disorder is to scatter the quasiparticles, 
and thus broaden any Fermi surface features (i.e., nested 
segments) which might give rise to peaks at particular 
q's in ^(q, w). Thus, in weak coupling, we expect weak 
quenched disorder to suppress, rather than enhance low 
frequency signatures of stripe order or incipient stripe or- 
der. Distinct responses to an applied magnetic field can 
also be expected to differentiate the weak and strong cou- 
pling limits. In particular, where the zero-field ground 
state is superconducting, induced order in vortex cores 
occurs under broad circumstances in the strong-coupling 
limit where there are competing orders; this sort of phe- 
nomena is generally much less prominent in the weak 
coupling limit. These issues are discussed in a number of 
good recent reviews - see footnotes 23 and 24. 



3. Intermediate coupling 

It is worth remembering that in many materials, the in- 
teraction strength is comparable to the Fermi energy, i.e., 
stability of matter (or some memory of the Virial theo- 
rem) conspires to place systems in the awkward regime of 
intermediate coupling. Here, neither the weak nor strong 
coupling approaches are well justified, and one is typi- 
cally forced to extrapolate results beyond their regime of 
validity. In this case, it is generally sensible to study the 
problem from both the strong and weak coupling per- 
spectives. Some features (such as those that are most 
sensitive to the presence of a particle-hole continuum) 
may be best viewed from the weak coupling perspective, 
while others (such as the role of collective modes) may 
be better viewed from strong coupling. 

One traditional way to understanding this regime is to 
include the effects of interactions in the context of the 
random phase approximation (RPA). For the complex 
spin susceptibility x, one writes 



X(k,w) 



Xo(k,w) 



1 - J7kXo(k,a;) 



(6.3) 



where xo is the susceptibility for noninteracting elec- 
trons, and U characterizes the interaction. In the case of 
dominantly forward-scattering interactions (such as the 
long-range piece of the Coulomb interaction), the RPA 
summation of bubble diagrams is justified even when the 
interactions are relatively strong. For the case of short- 
range interactions, no such systematic justification exists 
to the best of our knowledge. However, the results ob- 
tained from RPA are explicit and often intuitively appeal- 
ing. It is certainly reasonable, at least on a phenomeno- 
logical level, to compare the results of an RPA treat- 
ment with experiment to determine whether the physics 
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is simply connected to the weakly-interacting limit, or on 
the other hand whether the interactions produce quali- 
tatively new phenomena. 

One of the first RPA calculations was performed by 
Bulut et al. (1990) to model the antiferromagnetic spin 
fluctuations detected by NMR; for the interaction, they 
selected Uk = U, the on-site Coulomb repulsion. Lit- 
tlewood et al. (1993) found significant structure in the 
bare xo" (k, ui) from nearly-nested features of the Fermi 
surface. By suitable choice of the chemical potential, 
and inclusion of band narrowing by a factor of 2-4 com- 
pared to band structure calculations, they were able to 
find reasonable agreement with the early inelastic neu- 
tron scattering measurements of low-energy magnetic ex- 
citations in La2-£!;Sra;Cu04 (Chcong et al, 1991). Calcu- 
lations for YBa2Cu306+j, (Monthoux and Pines, 1993; Si 
et al., 1993) have typically found it necessary to employ 
the RPA expression with f/k = —J[cos{kxa) + cos{kya)] 
(or a similar form) in order to obtain reasonable con- 
sistency with experiment. Several recent calculations 
(Brinckmann and Lee, 1999; Norman, 2000; Onufrieva 
and Pfeuty, 2002) have addressed the resonance peak in 
YBa2Cu306_|_j, and the downward-dispersing excitations 
that appear at slightly lower energies (Aral et al., 1999; 
Bourges et al, 2000; Ito et al, 2002). 

In addition, we note that the lack of dispersion of the 
spin gap that develops below in La2-a;Sra;Cu04(Lake 
et al., 1999), together with the enhancement of intensity 
above the gap (Mason et al., 1996), have been interpreted 
in terms of the spin response of the electrons that par- 
ticipate in the superconducting state (Lake et al., 1999; 
Mason et al, 1996; Morr and Pines, 2000). 



B. Experimental evidence supporting a strong-coupling 
perspective in the cuprates 

From a purely theoretical perspective, it may not be 
possible to determine a priori whether there is an advan- 
tage to either a weak- or strong-coupling approach; how- 
ever, we believe that numerous experimental results point 
decisively towards the strong-coupling picture. Many of 
the relevant results involve measurements of magnetic 
correlations, and we have collected selected experimen- 
tal quantities in Table I. (Some explanations of how the 
parameter values are taken from experiment are given 
in App. B.) Below we highlight some of the important 
observations, and explain their significance. 

• Magnetism in undoped cuprates is quantitatively 
consistent with super-exchange between local mag- 
netic moments on copper ions. 



The undoped parent compounds, such as La2Cu04 
and YBa2Cu306 are antiferromagnetic insulators with a 
charge-transfer gap of 2 eV (Kastner et al, 1998). The 
ordered magnetic moments are consistent with one un- 
paired spin per planar Cu ion, after zero-point spin fluc- 
tuations, given accurately by spin- wave theory (Igarashi, 
1992), are taken into account. The effective super- 
exchange energy, Jgfj ~ 0.1 eV, determined from mea- 
surements of the spin-wave dispersion, has been calcu- 
lated from ab initio cluster models (Muhoz et al., 2000; 
Van Oosten et al., 1996), consistent with the local super- 
exchange mechanism. 

If the antiferromagnetic order corresponded to a spin- 
density wave due to Fermi-surface nesting, one would ex- 
pect the magnetic correlations and the optical gap to dis- 
appear at the Neel temperature. Neither of these things 
happens; these materials remain correlated insulators in 
the disordered state, and the spin correlations are well de- 
scribed by a nonlinear sigma model (Chakravarty et al., 
1989) with parameters derived from the ordered state 
(Kastner et al, 1998).^^ 

• With light doping, local magnetic moments change 
little, even though long-range Neel order is de- 
stroyed. 

In a /LtSR study of lightly-doped La2_2:Sr2;Cu04 and 
Yi_a;Caa;Ba2Cu306, Niedermayer et al. (1998) showed 
that there is essentially no change in the low-temperature 
ordered moment per planar Cu as the long-range antifer- 
romagnct order is destroyed by increasing the hole con- 
centration. This behavior has also been detected in NMR 
studies (Chou et al, 1993). In the "spin-glass" regime 
that occurs in the doping range between the antiferro- 
magnetic and superconducting phases, the ordered mo- 
ments decrease gradually (Niedermayer et al, 1998). We 
now know that for La2-i;Srj;Cu04 with 0.02 < x < 0.06 
the low-tcmperature phase involves ordering of magnetic 
moments in a diagonal spin-stripe structure (Fujita et al, 
2002; Wakimoto et al, 1999). That identical behav- 
ior occurs in Yi_a,Caa:Ba2Cu306 is supported by an 
NQR/NMR study (Singer and Imai, 2002). 

• The energy scale of magnetic excitations and the 

strength of the dynamic structure factor at frequen- 
cies less than J change only modestly with doping. 

One can see from Table I that ,/cfr decreases a rela- 
tively small amount with doping. There is also evidence 
in some underdoped YBa2Cu306+y samples for effects 
reminiscent of spin- wave- like dispersion at energies above 
50 meV (Bourges et al, 1997b; Fong et al, 2000). The 



Anderson (1997) has made a strong and compelling argument 
that antiferromagnetic exchange cannot be properly treated in 
RPA, and that the approach is internally inconsistent. 



^■"^ The recent observation by Coldea et al. (2001) that the spin-wave 
dispersion at large energies is better fit if a four-spin exchange 
interaction is included in the microscopic Hcisenberg model is 
still consistent with a system of localized spins. 
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TABLE I Various quantities characterizing the strength of antiferromagnetic correlations in doped and undoped cuprates and 
nickclatcs. All arc determined at low temperature. The quantity p is the nominal hole concentration per planar Cu; for 
YBa2Cu30fi+y and YBa2Cu408 this was estimated from the "universal" formula relating Tc/T^^^ and p from Presland et al. 
(1991). The acronyms used under "phase" are: AFI = antiferromagnetic insulator, CSO = charge-stripe order, SSO = spin 
stripe order, SC — superconducting order, rn is the magnetic moment per planar Cu, and m/mo is the moment relative to the 
AFI parent. J S(k, ui) is the dynamical structure factor, integrated over a Brillouin zone and over energy up to 100 meV, and 
normalized to the AFI parent. Jes is the effective super-exchange energy characterizing the maximum energy of the magnetic 
excitations (2Jeff ) measured by neutron scattering, or obtained from Raman scattering assuming that the two-magnon peak is 
at 2.7Jeff. Further discussion is given in Appendix B. 



material 


p 


phase 


m 
(Mb) 


m/mo 


rS(k,a;) 


Jeff 

(meV) 




^crit 


La2Cu04 





AFI 


0.60(5)" 


1 


1 


"> — n; 

152(6)" 




i 300 K 


La2Cuo.8(Zn,Mg)o.2Cu04 





AFI 




0.72(5)^^ 








i 100 K 


Lai.88Bao.i2Cu04 


0.12 


CSO, SSO 




0.60(2)"* 






^SSO 


Rj 35 K 


Lai.48Nd().4Sro.i2Cu04 


0.12 


CSO, SSO 




0.60(2)'* 




120(25)^^ 




« 50 K 


La2CuO4.11 


0.16 


SSO, SC 




0.62(2/ 








= Tc = 42 K 


Lai.86Sro.i4Cu04 


0.14 


SC 








0.43 


140(10)s 


Tc = 


35 K 


YBa2Cu306.i 





AFI 


0.55(3)'' 


1 


1 


125(8)* 




J 400 K 


YBa2Cu306.5 


0.09 


SC 








0.5-' 






52 K 


YBa2Cu306.6 


0.10 


SC 








0.6* 


110(15)* 


Tc = 


63 K 


YBa2Cu306.7 


0.11 


SC 








0.3' 


125(5)"" 


TcW 


70 K 


YBa2Cu408 


0.13 


SC 










125(5)"" 


Tc = 


80 K 


La2Ni04 





AFI 


1.0(1)" 


1 




31(1)° 


Tn R 


i 330 K 


La2Ni04.i33 


0.27 


CO 


0.8(1)'' 


0.80(5) 






Tsso 


= 110 K 


La2-a:Sra;Ni04 


~ 0.33 


CO 




1.11(5)^ 




20(1)'' 


Tsso 


« 200 K 



"Prom neutron diffraction (Yamada et al., 1987). 

"Prom inelastic neutron scattering by (Coldca et al., 2001). 

'^Prom neutron diffraction (Vajk et al., 2002). 

"^Local moment relative to La2Cu04 from ^SR (Nachumi et al., 

1998) . 

Prom two-magnon Raman scattering, assuming that the observed 
peak is at 2.7J (Gozar et al, 2003; Nachumi et al, 2002). 

^Local moment relative to La2Cu04 from /iSR (Savici et al, 
2002). 

sprom inelastic neutron scattering (Hayden et al., 1996a). 

''Prom neutron diffraction (Casalta et al, 1994). 

'Prom inelastic neutron scattering (Hayden et al, 1996b) and 

two-magnon Raman scattering (Blumbcrg et al, 1996, 1994). 

^Bascd on inelastic neutron scattering results (Fong et al, 2000; 
Hayden et al, 1996b). 

*Based on inelastic neutron scattering results (Dai et al, 1999; 
Hayden et al, 1996b). 

'Based on inelastic neutron scattering results (Fong et al, 2000; 
Hayden et al, 1996b). 

""From two-magnon Raman scattering (Blumberg et al, 1996, 
1994). 

"From neutron diffraction (Wang et al, 1992). 

"From inelastic neutron scattering (Yamada et al, 1991). 

pProm neutron diffraction (Tranquada et al, 1995a). 

^Local moment relative to La2Ni04 from /iSR (Jestadt et al, 

1999) . 

'^From inelastic neutron scattering (Boothroyd et al, 2003; 
Bourges et al, 2003). 



partially-integrated dynamical structure factor decreases 
with doping, but remains substantial. 

If the magnetic excitations corresponded to electron- 
hole excitations across the Fermi surface, one would ex- 
pect that the magnetic energy scale should increase to 



a value comparable to the Fermi energy. For that scale 
to be the same as JeS of the undoped phase would seem 
to be an incredible coincidence. With a substantial in- 
crease in energy scale, one would expect a corresponding 
decrease in the integrated dynamical structure factor. 
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Similar trends in terms of a modest reduction of Joff 
with doping and the survival of spin- wave- like excitations 
are found in stripe-ordered La2-a;Si'a;Ni04 (Boothroyd 
et al., 2003; Bourges et al., 2003). This system is def- 
initely in the strong-coupling limit, as it remains semi- 
conducting above the charge-ordering temperature (Kat- 
sufuji et al, 1996). 

• The k dependence of the magnetic neutron scat- 
tering measured from superconducting samples of 
YBa2Cu3 06+y is consistent with spin density on 
planar Cu ions. 

The magnetic scattering cross section is proportional 
to the square of the magnetic form factor (Fourier trans- 
form of the magnetization density); this is true for 
both elastic and inelastic scattering. An early analy- 
sis of inelastic magnetic scattering at different k points 
in superconducting YBa2Cu306+i/ samples revealed an 
anisotropy (Rossat-Mignod et al., 1992) that was later 
shown to be consistent with the anisotropy due to the 
magnetic form factor for a single spin in a 3dx2-y2 or- 
bital (Shamoto et al, 1993). 

The coherence of spin excitations within the Cu02 bi- 
layers of YBa2Cu306+j/ also results in a bilayer structure 
factor. The structure factor is a sinusoidal function of 
momentum transfer perpendicular to the bilayers, with 
a period that is proportional to the spacing between the 
layers of the centers of spin density. The observed period 
is quantitatively consistent with the intra bilayer spacing 
of the Cu atoms, but is incompatible with the spacing of 
the oxygen atoms, which is significantly smaller (Tran- 
quada et al., 1992). 

For electronic states at the Fermi level, there is a 
substantial amount of weight from planar O 2p states 
(Takahashi et al., 1988). As a result, one would ex- 
pect that magnetic scattering due to electron-hole ex- 
citations should, in real space, have significant weight 
associated with the oxygen sites. The measured modula- 
tions (Bourges et al, 1996; Fong et al., 1996; Tranquada 
et al., 1992) are not compatible with a significant com- 
ponent corresponding to the oxygen spacing. 

• When stripe order is observed in more heavily 
doped samples, the magnetic moments are large. 

As indicated in the table, the local magnetic moments 
detected in stripe-ordered phases by /xSR are comparable 
to the moment found in undoped La2Cu04. The large 
moments imply a strong modulation of hole density. The 
magnetic form factor determined by neutron diffraction 
is consistent with that expected for spin moments (Tran- 
quada et al., 1996). 

• When stripe order is observed, charge orders before 
the spins, and the ordering wave vector grows with 
doping. 

When charge and spin-stripe order are observed, 

as in Lai.6-a;Ndo.4Sr2;Cu04, the charge order ap- 
pears at a higher temperature than the stripe order 



(Ichikawa et al, 2000; Tranquada et al., 1996). (In 
La'i.875(Ba,Sr)o.i25Cu04, the ordering temperatures for 
charge and spin are very close (Fujita et al., 2002); nev- 
ertheless, the charge order parameter grows more rapidly 
than the spin order parameter.) It follows (Zachar et al., 
1998) that the charge order is not driven by the spin 
ordering. 

This result is natural for a strong-coupling picture of 
stripe correlations, but it poses a considerable challenge 
for the weak-coupling approach. In the latter, if one at- 
tributes the spin-stripe order to a nesting instability of 
the Fermi surface, then one must find a distinct nesting 
feature, with a spanning wave vector corresponding to 
Qchj to explain the charge-stripe order. Of course it is 
also necessary to explain the doping dependence of the 
ordering wave vectors. According to ARPES studies, the 
spanning wave vector near k = (27r/a)(|,0) is close to 
Qch- However, its variation with doping is opposite to 
that of Qch (Ding et al, 1997; Ino et al, 2002); with in- 
creasing X, the hole-like Fermi-surface approaches closer 
to the (tt, 0) point, so the spanning wave vector decreases, 
while Qch increases. We are not aware of a plausible ex- 
planation for the combined charge and spin modulations 
from the weak-coupling perspective. 

• Substituting Zn impurities into the Cu02 planes 
pins stripe order in La2_a;Sra;Cu04. 

The substitution of Zn into the Cu02 planes causes 
a significant increase in normal-state resistivity without 
modifying the carrier density (Fukuzumi et al, 1996). 
It also wipes out the dispersive feature near the nodal 
point observed by ARPES in the normal state (White 
et al, 1999). We have already discussed the fact that Zn- 
doping can induce static stripe order in La2-2;Sra;Cu04 
(see Fig. 8). These results arc incompatible with a mech- 
anism for stripe order based on Fermi-surface nesting. 
The Zn impurities break translational symmetry and 
cause considerable scattering of the charge carriers. Any 
sharp features at the Fermi surface are smeared out in 
the presence of the Zn. Thus the Zn should destroy, not 
induce, a nesting instability. 

• There is no evidence of well-defined quasiparticles 
in the normal state of under- and optimally-doped 
cuprates. 

There arc many features of the doped system which 
suggest that, for the most part, there are no well defined 
quasiparticles. ARPES spectra in Bi2Sr2CaCu208+5 and 
La2-2;Sra;Cu04 in some cases do exhibit features with a 
well defined dispersion, and by looking at the spectral in- 
tensity integrated over a small energy window about the 
Fermi energy, a Fermi surface of sorts is observed (Dessau 
et al, 1993; Ding et al, 1996; Ino et al, 2002). However, 
in no case has a peak been observed with a width small 
compared to its mean, which we take to be the definition 
of a well defined quasiparticle. In some ranges of temper- 
ature and k, there is no well defined quasiparticle peak, at 
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all. Near the nodal points, there are marginally defined 
quasiparticle, in the sense that there is certainly a clear 
peak in the spectral function, but its width is approxi- 
mately twice its mean energy (Valla et al., 1999). There 
is corroborating evidence that any quasiparticles are at 
best "marginal" which comes from the T-linear depen- 
dence of the normal state resistivity, and various other 
indirect (but bulk) measurements (Anderson, 1987, 1992; 
Lee, 1999; Varma et al, 1989). 

In addition, there is considerable indirect evidence that 
the familiar Fermi liquid power laws are strikingly absent. 
NMR experiments reveal significant temperature depen- 
dence to the nuclear 1 /TiT and dramatic violations of the 
Korringa relation, not only in the pseudogap regime, but 
even in the high temperature regime at optimal doping 
(Takigawa et al, 1991). The optical conductivity has a 
clearly non-Drude form (Collins et al., 1989), above and 
beyond the peculiar temperature dependence of the DC 
conductivity. The Hall number, as well, is anomalously 
temperature dependent (Wang et al., 1987). 

All together, these results strongly imply that a picture 
of weakly interacting quasiparticles has limited validity in 
the cuprates. The one exception to this concerns the low- 
energy behavior (i? <C Aq) of the cleanest superconduct- 
ing materials, deep in the superconducting state. Here, 
considerable indirect evidence exists that there arc re- 
markably long lived nodal quasiparticles which dominate 
the physics (Sutherland et al., 2003). However, direct 
evidence of these sharply defined excitations has yet to 
emerge in any single-particle experiment such as ARPES 
or STM. 



VII. CONCLUSIONS 

There are many reasons why identifying and study- 
ing fluctuations associated with the presence of "nearby" 
ordered states has become one of the main thrusts in 
the study of cuprate superconductors and related mate- 
rials. To some extent, these states are interesting just 
because they occur somewhere in the (multidimensional) 
phase diagram of these intriguing materials. Quantum 
critical points associated with some of these orders have 
been proposed to be the explanation of the notoriously 
peculiar high- temperature ( "normal state" ) behavior ob- 
served in many experiments. The possibility that certain 
local orders are inextricably linked with the phenomenon 
of high-temperature superconductivity is the most en- 
ticing reason of all. However, before we can determine 
whether a particular form of local order, such as local 
stripe order, could possibly be central to the problem of 
high temperature superconductivity, we need to deter- 
mine whether or not it is ubiquitous in the families of 
materials which exhibit high-temperature superconduc- 
tivity. Regardless of one's motivation, it is clear that the 
ability to identify the proximity to an ordered phase by 
the signatures of "fluctuating order" is very useful. This 
paper has focused on practical considerations pertinent 



to this task. 

A combination of rather general scaling considerations 
concerning quantum critical phenomena, and specific in- 
sights gleaned from the solvable models studied, has lead 
us to articulate a number of "lessons" concerning the op- 
timal way of obtaining information about nearby ordered 
states from experiment. 1) The information is best ob- 
tained from the low-frequency part of the dynamic struc- 
ture factor, preferably integrated over a small, but non- 
zero range of frequencies, with the scale of frequencies set 
by the characteristic frequency of quantum fluctuations, 
Ec/fi. 2) Weak disorder can make it easier, especially 
for static probes, to image the local order, as it can pin 
the fluctuations without greatly disrupting the intrin- 
sic correlations. 3) Experiments which reveal strongly 
dispersing features generally give information about the 
elementary excitations of the system; however, distin- 
guishing dispersing features that arise from well defined 
quasiparticles from the multiparticle continua character- 
istic of quantum critical points can be excecHlingly sub- 
tle. Specifically, even where one set of experiments can 
be sensibly interpreted in terms of band structure effects, 
care miist be taken in interpreting this as evidence of well 
defined electron-like quasiparticles. 4) Several aspects 
of the E and k dependence of the LDOS (measurable 
by STM) in the presence of weak disorder allow one, in 
principle, to distinguish interference effects due to the 
scattering of the elementary excitations from impurities, 
from the effects of pinned incipient order. Interference ef- 
fects in 2D produce peaks along curves in k space which 
disperse as a function of energy in a manner which is 
directly related to the quasiparticle dispersion relations, 
such as could be measured in ARPES, and they may (or 
may not) have a strongly energy dependent phase. Pin- 
ning of incipient order produces peaks at well defined 
points in k space which depend only weakly on energy, 
and generally have an energy independent phase. 5) It 
will often be true that interference effects and collective 
pinning will jointly produce complicated k and E depen- 
dent properties in the local density of states that arise 
from a combination of both effects, especially in energy 
ranges in which the interference and pinning features lie 
at nearby values of k. 

In Section IV of this paper, we applied these ideas to an 
analysis of the evidence of local stripe order in a number 
of neutron scattering and STM measurements on vari- 
ous cuprate superconductors. The evidence of both spin- 
and charge-stripe order is imambiguous in the La2Cu04 
family of high-temperature superconductors. However, 
there is increasingly strong evidence of substantial local 
charge-stripe order, and probably nematic order as well, 
in Bi2Sr2CaCu208+(5 and YBa2Cu306_|-y. Conversely, no 
evidence of incipient stripe order has been reported to 
date in the electron-doped cuprates, Ne2-xCexCu04±5 
and Pr2_xCea;Cu04±5. 

Note Added: After this work was completed, we re- 
ceived an advanced copy of a review article by (Sachdev, 
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2003) which has overlapping material with the present 
paper, and which reaches similar conclusions concerning 
the effects of incipient stripe order on the STM spectra. 
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APPENDIX A: Luttinger liquids as one-dimensional 
quantum critical charge ordered states 

In Section III we disciissed how a single impurity in- 
duces charge order in the simplest and best understood 
quantum critical system, the Tomonaga-Luttinger Model 
of a one-dimensional electron gas. In this Appendix we 
give a summary of the physics and of essential technical 
aspects of the discussion of Section III. In this Appendix 
we will only present the aspects of the theory relevant 
to the derivation of the expressions used in Section HI. 
There are a number of excellent reviews which cover the 
theory of the IDEG at great depth and we refer the in- 
terested reader to that literature (Emery, 1979; Fradkin, 
1991; Gogolin et al, 1998; Stone, 1994). 

Consider a one-dimensional system of interacting spin 
one-half fermions (electrons). We will denote by '^a{x) 
the Fermi field for an electron with spin a =t,i, and 
by ^±,cr(a;) its right and left moving components respec- 
tively: 

*,(ar) = e^^F^ V+,a(a;) + 6"^^^^ V-,<.(a;) (Al) 

where fcp is the Fermi wave vector. In what follows we 

will assume that the electron density is incommensurate 
and ignore umklapp scattering effects. Thus we will be 
working in a regime in which the dynamics of the right 
and left moving components of the electron are slowly 
varying and hence are well described by an effective con- 



tinuum Hamiltonian density H = Ho + Hmt where 
Ho = -ivF J2 {'>P+,M+<- - i^UM-,.) (A2) 

is the Hamiltonian density for non-interacting electrons 
in the low-energy regime where the dispersion is lin- 
earized; here v-p is the Fermi velocity. As visual, irrelevant 
operators which account for corrections to the linear dis- 
persion are not included (see however below for caveats). 
The effects of interactions are included in Jimt ■ 

The best way to describe the physics of the IDEG is 
by means of bosonization methods. Bosonization is the 
statement that the low-energy spectrum of the IDEG is 
exhausted by a long lived bosonic excitation described by 
a field 4>a-{x) which represents particle-hole fluctuations 
of spin a near the Fermi points ifep. It turns out that 
the field <pa also describes the phase fluctuations of a 2fcF 
CDW. Consequently the electron density operator Pa{x) 
is decomposed into a long wavelength piece jg and a 2A;f 
piece related to the density wave order parameter. The 
long wavelength electron density for spin projection <j, 
3o{x), is: 

Jo"G^)-<.^+,<.+V'l.V'-,a (A3) 

and the long wavelength current density with spin pro- 
jection cr, ii(;x), is: 

ilix) = - (A4) 

Here Jq and are operators normal-ordered with respect 
to the filled Fermi sea. 

The density and current density operators obey the 
equal-time bosonic commutation relations 

bo (^), jT(y)] = -d.5{x - y)5,,. (A5) 

TT 

where a, i/ =t, i. One can identify jo{x) and ji{x) with 

io (a;) = d^(t^a{x) 

3l{x) = U,{x) (A6) 

where (f)a{x) is a bose field and Ila{x) is its canonically 
conjugate momentum; they obey the canonical equal- 
time commutation relations 

[0^(x), H,.(y)] = iSa,^ 6ix - y). (A7) 

The (slowly varying) Fermi fields ip±,a{x) can themselves 
be reconstructed from the bose field (l)a{x): 

V±,<.(x) = Ma e^v^(±'/'-(^) - ^-(^)) (A8) 

(here a ~ v-p/D is a short-distance cutoff and D is the 
fermion band width) where we have introduced the dual 
field Ocr{x), defined by the identity 

dJa{x)=Iia{x) (A9) 
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and A/'cr, the Klein factor, is an operator which guarantees 
that Fermi fields with different spin labels anti-commute 
with each other. 

It is convenient to introduce the charge and spin bose 
fields Sr and 6^, 



(AlO) 



The Hamiltonian density of the interacting system can be 
written a sum of operators which are marginal (with scal- 
ing dimension equal to 2), relevant (with scaling dimen- 
sion smaller than 2) and irrelevant (with scaling dimen- 
sion larger than 2). For an incommensurate IDEG the 
effective low-energy Hamiltonian, which contains only 
marginal operators, has the spin-charge separated form 



W = Tic + 



(All) 



where the charge Hamiltonian density Tic in bosonized 
form has the universal (Tomonaga-Luttinger) form 

= ^K, + ^ {d,cp,f (A12) 

where Vc is the charge velocity and Kc is the charge Lut- 
tinger parameter. T For repulsive interactions, the spin 
Hamiltonian density Hg has the same form 



Hs = f H? + ^ (a.0.)^ 



(A13) 



The non-universal charge and spin velocities and the Lut- 
tinger parameters encode the dependence of this low en- 
ergy theory on the microscopic parameters of the sys- 
tem. Quite generally, for a system with repulsive interac- 
tions the charge Luttingcr parameter obeys the inequal- 
ity Kc < 1, and for spin rotationally invariant interac- 
tions the spin Luttinger parameter satisfies the equality 
Ks = 1. Typically the charge and spin velocities satisfy 
the inequality Vc > Vg- In the weak coupling limit, and 
neglecting all irrelevant operators, simple expressions for 
Kc, Vc and in terms of the backscattering and forward 
scattering amplitudes can be written down. 

Pliysic;al obscrvables of the IDEG have simple expres- 
sions in the bosonized theory. Since the Tomonaga- 
Luttinger model is strictly quadratic in bose fields, it 
allows for a straightforward computation of the correla- 
tion functions of all observables of interest, both at zero 
and finite temperature, as well as for different types of 
boundary conditions. For the purposes of this review it 



■'^ At intermediate and strong coupling the effective low-energy the- 
ory still has the same form but with significant renormalization 
of these parameters. 



will be sufficient to note that the electron density has the 
decomposition 

p{x) = +Mx) 

TT 

+ e^4fcF^O4fcp(a;)+e-^4fcFa.0t^^(^^ 



where 



30= Jo +jo 



IT 



(A15) 



(A16) 



is the long wavelength charge density, and 

= — cos(%/2^</.«)e-^^'^- 

is the CDW order parameter {i.e. , the 2fcF amplitude 

of the charge density), and 04kp is the Akp CDW order 
parameter (which we will not discuss here). 

In particular this implies that the charge dynamical 
structure factor for wave vectors close to 2kY, discussed 
in Section III, is the (retarded) correlation function of 
the CDW order parameter: 

ScDw{x,t) = {Ol,^^{x,t)OcDw{0,0)) (A17) 

whereas for small wave vectors it is given instead by the 
(retarded) density correlator 



So{x,t) = {jo{x,t)jo{0,0)) 



(A18) 



In particular, the spectral function for Scdw{x, t), which 
we will denote by ScDw{k,uj), has the scaling form 



ScDw{k,io) 



TTTa \ ^ / y h 



1677-^ V Vc 

(A19) 

where ^cDw{x,y) is the (dimensionless) scaling func- 
tion (Orgad et al, 2001) 



/OO POO 
du dv h 
-OO J — OO 



X h*Kc 



u — V 



hi - 



C30 

ru + V 



h\ 



U + V 

2 

ru — V 



(A20) 



Here r = Vg/vc < 1 and ha{z) is given by (Orgad et al., 
2001): 



hv{z) = Re 



{2iY B 



V ~ iz 



-,\-v 



where B{x,y) is the Euler Beta function: 

r{x)r{y) 



B{x,y) = 



Y{x + y) 



(A21) 



(A22) 
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and T{z) is the Gamma function. 

We now turn to the calculation of the tunneling density 
of states discussed in Section III. There wc presented the 
behavior of the 2kF component of the tunneling density 
of states N{2kF + q,E), the Fourier Transform in x and 
t, of the 2kF component of the electron spectral function. 
At zero temperature and for a semi-infinite system with 
< X < L (with L ^ oo), the 2kF component of the 
fermion Green's function, g2kF{xyt)j is given by (Eggert, 
2000; Eggert et al., 1996; Mattsson et ai, 1997) 

92kAx,t) = {il>l^^{x,t)^-^,{x,0)) (A23) 

^ 1 ^a+b-l/2 (^2x/VcTy 

27r {v^T - 2xf^'^ (vcT - 2xf {v^t + 2xf 



where t = t + iG+. Notice that in Eq. (A24) wc have 
kept the dependence on the short distance cutoff a; con- 
sequently, g2kj,{x^t) naively has units of a~^. However, 
in a Luttinger liquid the fermion operator has an anoma- 
lous scaling dimension. In for this semi-infinite geometry 
the scaling dimension is governed by the exponents a, h 
and c, which are given by 



{K, + 1)^ 



Sis:. 



At finite temperature T > 0, g2kF{x,t;T) becomes 



(A24) 



52fcp {x, t; T) 



1 f-KTa\ = /7rTa\ 
27ra \ Va J \ Vc J 



sinh 



(A25) 



^sinh (^^{vsT - 2a;)) J ^sinh (^(ucT - 2a;)) ] ^sinh (^^{vcT + 2a;)) j \ sinh (^Zr) 



r 



As with the structure factor, N{q + 2kF,E) can also be Kc and on the ratio of the charge and spin velocities 
expressed in terms of a scaling function $ as 



/•OO /-OO . . 

N{q + 2kF,E)= / dx / dt e-H92; - Et) g^^^ 

JQ J -OO 



Vc/Vg. B is the dimensionless quantity 
B = 



B faE 
2E 



2b 



T{a + b + c-^) 



(A27) 



2E E 



Vcq' ksT 



(A26) 



The behavior of iV(2fcF + q, E; T) for T > is shown 
in Fig. 3 and Fig. 4 in Section III. At T = the scaling 
which depends also on the charge Luttinger parameter function ^{2E/vcq) is given by 



$(u) 



Vs Jo 



dt t" (1 - t) 



a+b- 



{i-uty {i + u*ty 



1 - —Ut 



(A28) 



Here u = 2\E\/qVc + iO~^ ■ We have obtained analytic 
expressions for the singular pieces of N(k, E) ai T = 0, 
for a general value of the charge Luttinger parameter in 
the range of interest < Kc < 1, and general Vc/vs- 
N{k,E) has power law singularities as 2E/vcq ±1, 
2E/Vcq 0, and 2E/vcq oo and 2E/vcq Vg/vc- 

Physically the non-analyticities at 2E/vcq = ±1 and 
2E/vsq = 1 represent the threshold of the continuum of 
propagating excitations of appropriate type moving to 
the right and to the left respectively. The singularity at 
2E/vcq — > OO corresponds to the 2kF CDW stabilized 
(pinned) by the boundary. 

Close to the right-dispersing charge-related singularity 



at g — > 2E/vc we find 



1 faE 

Ni2kF+q,E)^-(- 



2b 



1 - 



2E 
Vcq 



(A29) 



where is a finite complex coefficient determined by 
the strength of the singularity. An important feature of 
this result is that as E goes through this singularity at 
fixed q, the phase of N{2kF + q, E) jumps by 7r(6— 1/2). 
(See Fig. 4.) Close to the right-dispersing spin-related 
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singularity at g ^ ^E/vs we find 
N{2kF + q,E) - 



by 



(A30) 



1 

2E 



oE 

Vs 



2b 



a+fa+c-l 

A' 



1 - 



2E 



2b- 



where A' is another complex coefficient. In the non- 
interacting limit Kc = 1 and Vs = Vc, these two singu- 
larities coalesce into a simple pole for a particle moving 
to the right. However, for Kc < 1, the TLL does not 
have quasiparticles but massless soliton states instead. 
As usual, the power laws reflect the multi-soliton contin- 
uum. (More generally, the limit Vg/vc — > 1 is somewhat 
subtle (Bindloss et al, 2003).) Similarly, for q -2E/vc 



N{2kF + q,E) 



1 /aE 



2E 



2b 

A^ 



2E 



(A31) 



where, once again. A- is a finite complex coefficient spe- 
cific to this singularity. However, unlike the singularity 
at E = +Vcq. the behavior near E = —Vcq although 
non-analytic is not divergent since & > for repulsive 
interactions. 

In addition to propagating excitations, we also find 
that N{2kp + q, E) has a singularity associated with the 
non-propagating CDW as g ^ 0; 



2b 



1/2 



2E 



Vcq 



(l--ffc)/2 



(A32) 

which diverges as g ^ 0. This singularity also exhibits a 
phase jump as g ^ 0''=, equal to 7r(l — Kc)/2. 
At low voltages £J — > and at fixed q, we find 



N{2kF + q,E)(x 



1 



aE\ 



2b 



2E \vc J 



1/2 r 



2E 



-I c+l 



Vcq 



(A33) 

Finally, we summarize here the calculation of the in- 
duced density of states in the weak-impurity limit, E > 
Tk- In Section II.C.l we showed that the density of states 
at wavevector k and voltage E induced by an impurity 
potential with amplitude V{k) at wave vector k defines 
the susceptibility XDOs{k, E) given in Eq. (2.8). Here we 
will sketch the calculation of XDOs(2fcF + q^E) (to lead- 
ing order in V{k)) for a Luttinger liquid. We must ffist 
find the Green's function 

G(z,i,Z2,23) = (A34) 
(V;,t(^i)V'-,t(^2) (V'i,T(^3)V'+,T(^3)+ T-i)) 

which can be computed readily using bosonization meth- 
ods. At T = and for an infinite system the time-ordered 
Green's function G{zi, Z2, z^), in imaginary time, is given 



G{zi,z2,Zi) = 



1 



\Zl - Z2\ 



Kc + 1 _1 

X {Zi - Zg) 4 (Z2 - Z3) " (Wi - W3) ^ 
X (Zl - 23)"^ {Z2 - Z3)~^ {W2 - W)3)"' 

(A35) 

where we have used the charge and spin complex coordi- 
nates z = X + iVcT and w = x + iVgT. The susceptibiUty 
XDOs{2kp + q, E) is obtained, upon analytic continuation 
to real frequencies, from 

dx / rfre-H'^a; - Et) 

-00 J —OO 

X j dr' G{x,t;x,0;0,t') (A36) 

The result for the integrated density of states N{2k-p + 
q,D) presented in Section III, Eq. (3.13) follows from 
integrating Eq. (A36) over energies large compared with 
the low-energy cutoff Tk- 



APPENDIX B: Experimentally determined scales of 
magnetism in various materials 

Here we explain some of the assumptions behind the 
parameter values specified in Table I. For antiferromag- 
netic phases, the meaning of the parameters is relatively 
straightforward, but extracting parameters from mea- 
surements on doped systems can require extrapolations 
from simple models. Note that we have not addressed 
important issues related to possible novel forms of or- 
bital magnetism that have been proposed (Chakravarty 
et al, 2001a; Varma, 1997). 



1. Magnitude of the ordered moments 

a. Absolute moments from neutron scattering 

In the case of long-range antifcrromagnctic order, the 
absolute magnitude of the average ordered moment per 
planar Cu, m, is obtained from an elastic neutron diffrac- 
tion measurement by comparing the intensities of the 
magnetic Bragg superlatticc peaks with those of Bragg 
peaks from the chemical lattice. It has been observed 
that m is correlated with the Neel temperature, In, in 
both La2Cu04 (Yamada et al., 1987) and YBa2Cu306+y 
(Tranquada et al., 1988). In the earliest measurements 
of La2Cu04, samples typically contained excess oxygen, 
resulting in values of m and Tn significantly lower than 
those of the stoichiometric material. Another complica- 
tion is that determination of m requires knowledge of 
the magnetic form factor, which turns out to be consid- 
erably more anisotropic (Shamoto et al., 1993) than was 
appreciated initially. 
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For comparison, it is worth noting that for a spin i 
Heiscnberg model on a square lattice with only nearest 
neighbor interactions, the magnitude of the average or- 
dered spin per Cu, to second order in a 1/S expansion 
(Igarashi, 1992), is (S) = 0.307. For a moment due to 
spin only, g = 2 and m = g{S) /ib = 0.61 ^b- Although 
it has not been possible to measure the g factor for Cu 
in the layered cuprates,^'' a typical value for a Cu^+ ion 
in a distorted octahedral environment is 2.2 (Abragam 
and Bleaney, 1970), which implies a typical moment of 
0.67 /iB- The fact that the observed maximum moments 
in La2Cu04 and YBa2Cu306+j/ are 10 to 20% smaller 
than this could be a result of hybridization of the Cu 
3dx2-y2 orbital with O 2pa orbitals. (By symmetry, the 
net spin density on oxygen sites is zero in the Neel struc- 
ture.) 

b. Relative moments from /xSR 

In /xSR, one measures a precession frequency, which 
is directly proportional to the local hyperfine field at 
the /i"*" site. Determination of the absolute moment re- 
quires knowledge of the location (obtained by cal- 
culation rather than experiment) and a calculation of 
the magnetic field at that site due to the ordered Cu 
moments. Alternatively, under the assumption that the 
/x+ site is unchanged by doping, one can determine the 
local moment in a doped sample relative to that in 
the parent insulator. An advantage of the /zSR tech- 
nique is that its sensitivity docs not depend on the exis- 
tence of long-range order. A beautiful example of this 
is the study of Niedermayer et al. (1998) in which it 
was shown that the local ordered moment (at very low 
tempcratiirc) changes little as long-range antiferromag- 
netic order is destroyed by doping in La2-xSr3;Cu04 
and Yi_a;CaxBa2Cu306. In a complementary study, 
Klauss et al. (2000) have shown that substantial moments 
associated with local antiferromagnetic order (presum- 
ably stripe order pinned by the tetragonal structure) are 
present in Lai.8-xEuo.2Sra;Cu04 for 0.08 < a; < 0.18. 

In a stripe-ordered sample, the muons sample a distri- 
bution of local hyperfine fields, resulting in a fairly rapid 
damping of the precession signal (Klauss et al., 2000; 
Nachumi et al, 1998; Savici et al, 2002). The relative 
moment quoted in the table corresponds to the maximum 
of that distribution. 



2. Integrated low-energy spectral weight 

Neutron scattering measures the dynamical magnetic 
structure factor S'(Q,a;) = S'zz(Q,a;). Actually, the 



Despite numerous attempts, it has not been possible to detect 
an electron-spin resonance signal from Cu^"*" in planar cuprates 
(Simon et al, 1993). 



spin structure factor is the tensorial quantity /S'a/3(Q, w), 
which obeys the sum rule 

[ dw f (iQVS„„(Q,w) = 5(5 + 1), (Bl) 
Jo Jbz ^ 

where S is the spin per magnetic ion. For cuprates in 
general, and superconducting samples in particular, it is 
challenging to measure the magnetic scattering to suffi- 
ciently high energies (with sufficient signal-to- noise ratio) 
that the sum rule can be evaluated; however, there do 
exist a few studies that allow one to evaluate the energy 
integral up to 100 meV. While there arc significant un- 
certainties in evaluating the integral from published data, 
these should be smaller than the uncertainty in the cali- 
bration of the sample volume contributing to the signal, 
which for superconducting samples can be on the order 
30% (Fong et al, 2000). 

Since we are interested in the change of this integrated 
quantity with doping, we have normalized it to the value 
measured for the relevant parent antiferromagnet. Note 
that for the latter case, the integral includes the signal in 
the magnetic Bragg peaks as well as that from the spin 
waves. 

Another quantity that one can evaluate is the relative 
weight under the resonance peak, which appears below 
Tc in the superconducting state. Values, normalized to 
the full sum rule weight, have been tabulated by Kee 
et al. (2002). For YBa2Cu306+y samples with vary- 
ing from 52 to 93 K, they found the normalized weight to 
be about 1.5%; for a sample of Bi2Sr2CaCu208-(-5 with 
Tc = 91 K, the value is ~ 6%. It should be noted that 
integrating S'(Q,a;) up to 100 meV for the experimen- 
tal measurements on antiferromagnetic YBa2Cu306.i5 
(Hayden et al, 1996b) gives just 20% of the full sum 
rule weight; hence, using the same normalization as for 
/ 5(Q, uj) in Table I, the weight of the resonance peak in 
YBa2Cu306-i-3/ is about 0.08. 



3. Energy scale of spin fluctuations 

a. From neutron scattering 

In the spin-wave theory of a Heisenberg model on 
a square lattice with only nearest-neighbor coupling, 
the zone-edge magnon with qo = (tt, 0) has energy 
= 45 Jeff = 2 Jeff, where the final equality applies for 
5 = i . Longer-range couplings have observable effects in 
the spin- wave dispersion curves measured for La2Cu04 
(Coldea et al, 2001); nevertheless, it is useful to char- 
acterize Eqo, which is essentially the local energy to flip 
a spin, in terms of Jcs- Assuming that a local super- 
exchange coupling is still active between local moments 
in doped samples, we have used the same formula to de- 
fine Jeff in terms of the maximum spin-excitation energy 
observed in superconducting cuprates. 
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b. From two-magnon Raman scattering 

While Raman scattering is not sensitive to individ- 
ual spin waves, it docs provide a valuable probe of two- 
magnon correlations in 2D square-lattice antiferromag- 
nets (Lyons et al, 1988). The dominant part of the re- 
sponse has Big symmetry and appears at low tempera- 
ture as a strong peak in intensity, with the peak occurring 
at an energy of 2. 8 Jeff for 5 = ^ and 6.8Jeft for 5 = 1 
(Canali and Girvin, 1992). The scattering mechanism is 
believed to involve relatively short-range excitations, so 
that the response is not strongly sensitive to the existence 
of long-range order. Although there is no formal justifi- 
cation for it, we have used the formula for peak energy in 
the antiferromagnetic state in order to extract a charac- 
teristic value for Jeff from the 2-magnon signal measured 
from superconducting samples. Note that there are also 
results available on Bi2Sr2CaCu208+5 (Blumberg et al., 
1997; Sugai and Hosokawa, 2000) which we have not in- 
cluded in Table I. 
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